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SECTION  1 
INTRODUCTION 

The  program  reported  here  was  initiated  to  improve  computer  modeling 
of  the  HCI-H2  non-reacting  mixing  laser  (NRML).  Figure  1 illustrates  the 
basic  laser  modeling  approach.  In  the  present  case,  the  equilibrium  chem- 
istry modeling  represents  the  combustion  processes  occuring  in  the  plenum 
where  a solid  fuel  is  burned  to  produce  HCl  and  H2.  Fluid  dynamics  rou- 
tines describe  the  flow  of  the  gases  through  the  throat  and  the  expansion 
through  a nozzle  into  the  laser  cavity  region.  Downstream  of  the  nozzle, 
non -equilibrium  chemistry  models  represent  the  HCl  pumping  reactions  and 
the  HCl  and  H2  deactivation  processes. 

The  resulting  specie  populations  are  then  used  with  the  basic  radiative 
properties  (molecular  linewidths  and  Einstein  coefficients)  to  determine  the 
incremental  gain  profile  of  the  medium.  Finally,  a cavity  model  is  em- 
ployed to  determine  the  power  and  spectral  content  of  the  laser  output. 

The  complexities  of  laser  device  models  vary  greatly  depending  on  the 
purpose  and  the  accuracy  required  in  a specific  application.  We  have  con- 
centrated on  developing  general  non-equilibrium  chemistry  and  radiative 
property  programs  which  can  be  basic  modules  in  either  simple  or  sophis- 
ticated laser  modeling  codes;  and  the  generation  or  compilation  of  specific 
input  data  for  modeling  the  HCI-H2  NRML. 

The  following  section  of  this  report  describes  the  linewidth  calcula- 
tions performed  for  this  program  and  a computer  subroutine  incorporating 
these  results  which  may  be  easily  included  in  most  laser  performance 
codes.  The  third  section  summarizes  the  current  status  of  Einstein  coeffi- 
cient modeling  for  HCl  lasing  transitions.  Next,  we  discuss  the 
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MERADCOM  non -equilibrium  chemistry  code  which  was  revised  and  ex- 
panded under  this  contract.  Our  summary  and  recc  mmendations  appear  in 
the  final  section  of  this  report.  The  appendices  contain  tabular  linewidth 
results  and  listings  of  the  codes  generated. 
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SECTION  2 

LINEWIDTH  MODELING 


The  basic  quantity  required  for  analyzing  a laser  system  is  the  local 
gain  of  the  medium.  This  is  given  as: 

hN 

a = -^w(PPB(2J  + 1) 


n(v  + 1) 


hcE 


exp 


v+1, J-1 
kT 


n(v)  


(1) 


The  notation  is  as  given  by  Emmanual  [1].  All  quantities  in  this  equation 
are  simple  functions  of  the  energy  levels  of  the  molecule  except  vibrational 
state  concentrations  (n(v),  n(v  + 1)),  the  line  profile  (<p),  and  the  Einstein 
coefficient  (B).  Of  interest  here  is  the  line  profile;  the  Einstein  coefficient 
is  discussed  in  a later  section,  and  the  state  concentrations  are  determined 
by  the  chemistry. 


The  detailed  line  profile  depends  in  a complicated  manner  on  the 
Doppler  and  Lorentz  linewidths  of  the  particular  laser  transition;  however, 
for  the  case  in  which  the  Lorentz  width  is  significantly  greater  than  the 
Doppler  width,  the  line  center  gain  value  is  simply  inversely  proportional 
to  the  linewidth. 


tr> 


LR 


(2) 


The  Lorentz  (or  collision-broadened)  linewidths  are  related  to  the  interac- 
tions of  the  lasing  molecules  with  other  molecules  in  the  medium.  Semi- 
classically,  Lorentz  broadening  results  from  collisions  undergone  by  the 


lasing  molecule  which  interrupt  the  lasing  process.  These  collisions  may 
be  with  either  similar  or  different  species  than  lasing  species.  In  the  case 
of  HCl  lasing,  other  HCl  molecules  are  very  effective  in  causing  Lorentz 
broadening  because  of  the  relatively  large  value  of  the  HCl  dipole  moment  and 
the  long-range  nature  of  the  diix>le -dipole  force.  On  the  other  hand,  homo- 
nuclear  molecules  such  as  H2  and  atoms  are  less  effective  broadening 
agents  because  they  have  no  dipole  moment. 

The  basic  equation  for  the  Lorentz  linewidth  is 


r(n>)  = (“'“j  (”>  (3> 

Jp  P P 

In  this  equation,  y is  the  half -width  for  a particular  rotational  transition 

designated  by  m.  The  critical  variables  in  (3)  are  the  rotational  populations 

of  the  perturbing  molecular  states  Pt  , and  the  broadening  cross-sections 

P 

of  the  perturbing  rotational  states  o t . For  the  present  discussion,  we  will 

P 

consider  vibrational  states  of  the  perturbing  molecules  as  separate  species. 


The  rotational  population  Pj  , is  only  a function  of  the  rotational 
energy  levels,  the  kinetic  temperature,  and  the  degeneracies  of  the  rota- 
tional states  of  the  perturbing  molecules.  The  critical  and  difficult  factor 
is  the  effective  broadening  cross-section,  Oj^.  This  can  be  expressed  as  [2] 


a = TT 
J 

P 


+ C,b“^2,g,.F,(k  ) + C„b’y  ,g„.F„(k  ) 
o 1 0 “*^11  1'  o 2 o — ‘^2i  2'  o 
I i 


^3^0  Z^^3i^3^o^ 

i 


(4) 


In  (4),  b^  is  an  impact  parameter  obtained  from  the  solution  of  the  follow- 
ing equation. 


(5) 


+ Cjb  ^ ' 

In  Equations  (4)  and  (5),  the  C.  are  coefficients  related  to  the  molecular 
moments  of  the  molecules,  g^.  represents  Clebsch-Gordan  coefficients  be- 
tween molecular  states,  and  the  are  resonance  functions.  The  F.  de- 
scribe strengths  of  various  interactions  as  functions  of  the  energy  defect  in 
the  rbtational  energy  exchange  in  the  collision  process.  The  F^  approach  1 
when  the  energy  defect  is  small  because  of  resonance  between  rotational 
transition  energies  in  the  emitting  and  perturbing  molecule. 


The  multipole  moment  parameters  are  contained  in  the  expressions. 


(6) 


(7) 

(8) 


where  /i.  and  Q.  are  the  dipole  and  quadrupole  moments  of  the  molecules, 
and  subscripts  1 and  2 refer  to  the  active  and  perturbing  molecules 
respectively. 


2.  1 HCl  SELF-BROADENED  LINEWIDTHS 

Values  for  Pj^  and  Cj^  calculated  for  the  self -broadening  of  the  HCl 
P2(5)  laser  transition  are  shown  in  Figure  2.  The  Uj  's  presented  as  bars 
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are  typical  of  self-broadening  where,  because  of  the  close  match  between 
the  rotational  energy  levels  of  the  molecules,  large  broadening  cross- 
sections  occur  for  J^'s  near  transition  J values.  Large  linewidth  values  re- 
sult when  large  cross-sections  coincide  with  large  values  of  the  distribution 
function.  For  moderate  J transitions,  the  self -broadened  cross-sections 
are  generally  much  larger  than  the  kinetic  theory  (or  hard  sphere)  cross- 
sections  of  the  molecules.  Compare  the  HCl  kinetic  cross-section  of 
~3  X 10  cm  with  the  broadening  cross-sections  plotted  in  Figure  2. 

Under  these  conditions,  the  Anderson  theory  used  here  is  reliable-  Fig- 
ure 3 shows  good  agreement  for  all  temperatures  measured  [3]. 

Since  the  initial  states  of  most  lasing  transitions  are  sparsely  popu- 
lated at  equilibrium  conditions  for  moderate  temperatures,  direct  measure- 
ment data  of  the  line -broadening  coefficients  for  these  transitions  are  diffi- 
cult to  obtain.  Thus,  knowledge  of  these  linewidths  must  be  based  on 
calculations.  A number  of  interesting  factors  come  into  play  when  we  con- 
sider these  linewidths. 


Temperature  dependence  of  individual  linewidths  is  influenced  mainly 
by  the  broadening  molecule's  population  distribution.  Low  temperatures 
result  in  a population  distribution  which  is  sharply  peaked  near  low  val- 
ues. For  the  self -broadening  cases,  this  means  that  large  Pj^  values  coin- 
cide with  large  values  at  small  m.  A distribution  of  linewidths  peaked 
for  small  m results,  as  in  the  curves  shown  in  Figure  3.  At  higher  tem- 
peratures, the  population  distribution  is  broader,  and  the  linewidths  have 
a weaker  m dependence. 


Since  the  rotational  constant  (B^)  of  a molecule  changes  as  a function 
of  the  vibrational  state,  the  rotational  resonances  shift  in  the  higher  vibra- 
tional states.  The  effects  of  these  shifting  resonances  on  the  linewidths 
are  shown  in  Figures  4 and  5.  The  widths  for  the  1 ->  0 bands  closely  follow 
the  shape  of  the  Boltzmann  distribution  function  pj^.  This  is  because  there  is 
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2.  CALCUi  ATED  POPULATION  Jp  UlSTHIBUTiON 
BROADENING  CflOSS- SECT  IONS  FOR  HCl  P2(5) 
BROADENED  BY  HCl  (v  - 0) 


FIGURE 

AND 


(wj,v/^_WD)  HxcnMaNn 


I 


9 


twite. ■., ,->:* a,- ?aMC.-~j<^.-->;A-v-...  »»«t.  j 

-5  -9  -13  -17  -2 

m 

FlGUliE  4 CALCULATED  LINEWIDTHS  FOR  THE  IICl 
(v  = 2 1)  LINES  BROADENED  BY  HCl  (v  = 0) 

AT  VARIOUS  TEMPERATURES 
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a close  match  between  the  resonances  in  the  broadening  molecule 
(HCl  (v  = 0))  and  the  active  transitions  (HCl  (v  - 1 -»  0)).  On  the  other  hand, 
a beat  phenomena  is  observed  in  the  linewidths  for  the  HCl  (v  = 2 -♦  1)  band 
because  of  the  shifting  resonances  between  these  transitions  and  the  broad- 
ening molecules  for  various  rotational  states. 

Similar  resonance  effects  occur  when  the  broadening  by  different  vi- 
brational states  of  HCl  are  considered.  This  can  be  seen  in  Figure  6, 
where  the  broadening  of  the  HCl  (2  1)  band  by  various  HCl  vibrational 

states  is  shown.  This  figure  illustrates  the  strength  of  resonance  interac- 
tion effects  for  low  m values:  factors  of  two  differences  are  calculated  for 
the  effectiveness  of  various  HCl  vibrational  states  in  broadening  the  same 
lasing  transition. 

2.  2 FOREIGN-BRO.ADENED  HCl  LINEWIDTHS 

The  principle  foreign  broadener  in  the  NRML  system  is  H^,  We  used 
the  basic  approach  described  in  2.  1 to  calculate  the  H^  broadening  of  HCl: 
however,  the  problem  is  not  as  straightforw-ard  as  with  self-broadening. 

The  difficulties  v/ith  a first-principle  calculation  of  these  foreign-broadened 
widths  are  that  the  interactions  are  weaker,  resonances  are  less  likely  to 
occur,  and  thus  the  broadening  cross-sections  are  comparable  to  or  smaller 
than  the  kmeiic  theory  cross-sections.  These  problems  cause  extreme  dif- 
ficulties with  the  Anderson  theory  because  the  small  broadening  cross- 
sections  imply  that  the  broadening  is  caused  by  collisions  with  impact 
parameters  less  than  the  molecular  diameters.  The  theory,  however,  as- 
sumes straight  line  trajectories  which  are  unrealistic  for  smiall  impact 
parameter  collisions.  An  alternate  method  for  treating  these  close  colli- 
sions is  critical  for  calculating  foreign  broadening  linewidths. 

Three  approaches  for  treating  these  close  collisions  are  illustrated  in 
Figure  7,  In  that  figure,  b is  the  impact  paraineter  of  the  collision,  and 
S(b)  is  the  collision's  effectiveness  in  interrupting  the  radiative  process. 
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FIGURE  6 CALCULATED  LINEWIDTHS  FOR  THE 
HCl  (v  = 2 - 1)  LINES  BROADENED  BY  VARIOUS 
HCl  VIBRATIONAL  SPECIES  AT  300  K 
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Anderson's  basic  approach  [4]  assumes  that  any  collision  with  an  impact 
parameter  less  than  the  kinetic  theory  diameter  is  totally  effective  in  caus- 
ing broadening.  Another  approach  using  the  Anderson  theory  is  simply  to 
ignore  the  unrealistic  straight  trajectory  approximation  for  close  collisions 
and  let  the  impact  parameter  he  as  small  as  necessary  to  satisfy  Equation  5, 
In  the  Van  Kranendonk  approximation  [5],  the  collision  efficiency  is  calcu- 
lated at  the  kinetic  theory  impact  parameter.  This  collision  efficiency  is 
then  assumed  for  all  impact  parameters  smaller  than  the  kinetic  theory 
value. 

HCl  linewidth  computations  based  on  these  three  assumptions  are  com- 
pared with  measurement  results  in  Figure  8.  As  illustrated  here,  the  basic 
Anderson  theory  predicts  larger  widths  than  were  measured.  The  Van 
Kranendonk  approach  shows  re”"  jonable  agreement  with  the  measurements 
for  small  and  large  m values  but  underpredicts  the  widths  for  moderate 
m's.  The  Van  Kranendonk  approach  fails  by  underpredicting  the  cross- 
sections  for  small  impact  parameter  collisions. 

To  fit  the  data,  we  use  the  Van  Kranendonk  approach  but  assume  that 
collisions  at  less  than  the  kinetic  theory  impact  parameter  are  at  least  40^c 
efficient  in  causing  broadening. 

The  linewidths  calculated  for  HCl  broadened  by  H2  with  the  above  as- 
sumptions are  shown  in  Figure  9,  Good  agreement  is  now  seen  between  the 
calculations  and  measurements.  The  predicted  hump  in  the  low  temperature 
linewidths  near  an  m of  -15  results  from  a resonance  between  the  radiation- 
less quadrupole  transition  in  K2  and  the  dipole  transition  in  HCl.  It  would 
be  an  interesting  test  of  the  present  theory  U experimental  data  could  be  ob- 
tained at  large  m value  where  the  hump  is  predicted.  Although  small  varia- 
tions in  the  linewidths  may  be  expected  for  different  vibrational  states  of 
the  active  and  perturbing  molecules,  we  have  not  calculated  them  in  the 
present  study. 
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FIGURE  8.  COMPARISONS  OF  EXPERIMENTAL  HCl 
BROADENED  LINEWEDTHS  WITH  VALUES 
CALCULATED  BY  VARIOUS  METHODS 
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FIGURE  9.  HCl  (1  -♦  0)  LINEWIDTHS  BROADENING 
BY  H2  (v  = 0)  AT  VARIOUS  TEMPERATURES 
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The  linewidths  calculated  for  300  K for  both  self -broadening  and  H2 
broadening  of  HCl  are  given  in  Appendix  L 


2.  3 LINEWIDTH  SUBROUTINE  DEVELOPMENT 

To  facilitate  the  inclusion  of  these  detailed  linewidth  calculations  in 
laser  modeling  codes,  a least -square  fitting  was  accomplished  of  HCl  line- 
widths  and  a subroutine  written  to  utilize  the  resulting  coefficients  in  radia- 
tive transfer  codes. 


The  series  of  exponential  functional  forms  listed  below  was  used  to 
represent  the  m dependence  of  the  linewidths. 

. V -m/4  -m/4  2 -m/2 

y(m)  = Cj  + C2  e + c^m  e + c^m  e 

-m^/8  2 -m^/16  -m^/8 

+ Cgm  e + Cgm  e + e 


(9) 


This  type  of  expansion  works  very  well  for  linewidths  [6].  Figures  10 
through  12  represent  the  quality  of  computer  fits  obtained  in  the  present 
study.  Figure  10  is  representative  of  the  better  fits,  Figure  11  is  a nom- 
inal match  and  Figure  12  is  one  of  the  poorer  fits.  The  poorer  ones  miss 
some  of  the  fine  detail  in  the  calculations  but  overall,  the  results  are 
reasonable. 


In  addition  to  fitting  the  HCl  self-broadening  and  H2  broadening,  ex- 
perimental linewidths  of  HCl  broadened  by  Ar  were  also  fit.  The  experi- 
mental widths  measured  by  Petrov  [7]  were  used  because  they  agree  with 
other  authors'  results  at  300  K [8]  and  also  include  higher  temperature  re- 
sults. We  plotted  the  widths  measured  by  Petrov  at  300  and  600  K on  log- 
log  paper  and  assumed  a straight  line  approximation  to  interpolate  or  ex- 
trapolate values  at  200,  400  and  500  K.  These  linewidths  were  then  fit  in 
the  same  manner  discussed  above  and  the  results  are  shown  in  Figure  13. 
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FIGURE  12.  LEAST-SQUA 
BROADENED  B’ 
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FIGURE  13.  LEAST-SQUARE  FITS  TO  HCl 
LINEWIDTHS  BROADENED  BY  Ar 


The  least-square  fit  parameters  for  all  HCl  widths  form  Appendix  IL 


2.  4 LINEWIDTH  EVALUATION  SUBROUTINE 

To  make  the  linewidth  results  easily  assessable  to  HCl  laser  modeling 
groups,  we  wrote  a computer  subroutine  which  uses  the  least-square  fit 
parameters  with  the  standard  chemical  variables  and  returns  the  computed 
pressure-broadened  linewidth.  Table  1 lists  the  variables  required  for 
calling  this  subroutine  and  gives  the  common  storage  assignments  for  the 
least -square  parameters.  The  basic  computation  done  by  the  program  is 
given  as  Equation  10. 

r(T,  m).p(  <*' 

i = perturbing  j = fitting 
molecular  states  functions 

Bookkeeping  complications  arise  because  the  laser  modeling  codes  may 
include  more  or  fewer  molecular  states  than  are  represented  in  the  least - 
square  coefficients,  and  temperature  interpolation  is  usually  required; 
however,  the  program  includes  the  more  general  cases.  Appendix  III  lists 
the  subroutine  GAMMA2  along  with  a calling  program  and  a simple  test 
case. 
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TABLE  1 


CALLING  SEQUENCE  AND  COMMON  ASSIGNMENT 
FOR  LINEWEDTH  SUBROUTINE  GAMMA2 

Calling  Sequence 

WIDTH  = GAMMA2  (PTOT,  TKEN,  CLS,  NCLS,  CNLS,  NCNLS,  CSA,  NCSA, 
LV,  M) 


Variable 

Description 

PTOT 

Total  Pressure 

TKEN 

Kinetic  Temperature 

CLS(  ) 

Array  of  the  mole  functions  of  the  vibrational  states  of  the 
lasing  species 

NCLS 

Number  of  mole  fraction  in  CLS 

CNLS(  ) 

Array  of  mole  fractions  of  the  vibrational  states  of  the 
non-lasing  species 

WCNLS 

Number  of  mole  fractions  in  CNLS 

CSA(  ) 

Array  of  mole  fractions  of  other  species 

NCSA 

Number  of  mole  fractions  in  CSA 

LV 

Vibrational  designation  of  lasing  transition 

M 

Rotational  designation  of  lasing  transition  (-Jfinal) 

Common 

Variable 

Description 

NPL 

Number  of  self -broadening  states  included 

NVL 

Number  of  vibrational  lasing  transitions  included 

NT 

Number  of  temperatures  included 

TS(  ) 

Temperatures  included 

COFT(  ) 

Self -broadening  least-square  fit  coefficients 

NPNL 

Number  of  foreign  broadening  vibrational  states 

CNOFT(  ) 

Foreign  broadening  coefficients 

NPA 

Number  of  other  broadening  species 

CAOFT(  ) 

Coefficients  for  other  species 
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SECTION  3 

HCl  EINSTEIN  COEFFICIENTS 


The  Einstein  coefficients  of  a molecular  transition  are  related  to  the 
electric  dipole  matrix  element  of  that  transition  by  the  following  formulae: 

4 3 

A(v' J'  ->  vJ)  = 3h(2.r  i”n  |p(r)  I vJ>  I ^ (moIecuIe-sec)'^ 


4 

B(v’J'  - vJ)  = ^ - l<v’J’  |p(r)|vj>l^ 

3h'^c(2j’  + 1) 

B(V  J'  vJ)  = ^^|j^-^^^B(v’J’  - vJ) 


cm 


.molecule -erg-sec 


Matrix  element  evaluations  for  laser  transitions  often  must  be  extrapolated 
from  matrix  elements  of  other  non-lasing  transitions  which  are  more  ac- 
cessible to  measurement.  This  procedure  was  followed  in  an  earlier 
series  of  calculations  by  the  present  authors  on  the  HCl  molecule  [9,  10]. 

In  this  earlier  work,  the  available  experimental  determinations  of  the  vi- 
brational matrix  elements,  and<0l^l2>,  were  reviewed  and 

"best"  values  selected. 


Since  publication  of  the  earlier  work,  we  have  become  aware  of  some 
additional  vibrational  matrix  element  values  obtained  from  dispersion  and 
absorption  measurements  of  gaseous  HCl  [11,  12,  13].  These  additional 
results  are  presented  in  Figure  14  along  with  the  measurement  values  in- 
cluded in  the  earlier  paper.  In  that  figure,  the  results  of  an  experimenter 
who  reported  values  of  both  matrix  elements  are  represented  as  a point. 
Where  a particular  author  has  measured  only  one  matrix  element,  this  re- 
sult is  presented  as  a slash  on  the  appropriate  axis. 
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FIGURE  14.  EXPERIMENTAL  DETERMINATIONS  OF  <0l>ill> 

AND<01m12>  for  HCl 
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The  results  of  Kimel  and  Legay  are  seen  to  agree  to  within  three  per- 
cent with  the  values  selected  for  use  in  the  earlier  calculations;  however, 
the  value  obtained  by  Boulet  is  nearly  6%  larger  than  the  value  of  ^0l/i|2> 
previously  assumed.  When  the  additional  data  is  considered  with  the  re- 
sults included  previously,  no  significant  trends  emerge.  Thus,  we  have 
decided  not  to  redo  the  Einstein  coefficient  calculations  in  the  present 
study. 


SECTION  4 

CHEMISTRY  CODE  DEVELOPMENT 


An  improved  set  of  finite  rate  chemistry  codes  was  developed  in  this 
program  from  previous  MERADCOM  codes.  A flow  chart  of  these  codes  is 
shown  in  Figure  15.  A list  of  important  variables  used  in  the  chemistry 
routines  is  given  in  Table  2.  As  presently  configured,  the  codes  treat  up 
to  99  species  including  catalytic  species,  and  accept  chemical  reactions  in 
a manner  very  similar  to  standard  chemical  notation.  The  codes  return 
the  enthalpy  change  from  the  chemical  reaction  and  the  rates  of  production 
of  the  various  species. 

The  calling  program  for  the  chemistry  routines  is  CHEMS,  CHEMS 
first  reads  the  thermodynamic  data  for  each  species.  These  inputs  are 
molecular  weight,  enthalpy  of  formation  at  25  K,  and  specific  heat  and 
enthalpy  as  a function  of  temperature.  Next,  the  reactions  to  be  included 
are  entered  and  followed  by  the  reaction  rate  coefficients.  A numeric  code 
designates  the  form  of  the  reaction  rates  to  be  used  for  each  reaction. 
Finally,  the  flow  variables,  pressures,  temperatures  and  time  increments 
are  entered,  and  CHEMS  calls  on  additional  subroutines  to  perform  the 
computations. 

The  housekeeping  routines  CHEMVT  and  CHEMSUM  compute  the  vari- 
ables required  for  the  V-T  deactivation  processes  and  for  using  the  com- 
posite compounds.  Then,  the  routine  CHMSTY  is  called  to  compute  the 
species  formation  rates. 

The  subroutine  CHMSTY  first  calls  the  MATH  subroutine  to  compute 
the  reaction  rate  at  the  specific  temperature  of  interest.  Next,  CHMSTY 
evaluates  the  individual  production  rate  for  each  species  and  each  reaction. 
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TABLE  2 

DEFINITIONS  OF  PARAMETERS  USED  IN 
CHEMISTRY  ROUTINES 


Variable 


Description 


AID 

ALPHA 

AT 

AVAGAD 

A,  Al,  A2 

AR 

CN 

CNU 

CNTI 

CPTB 

CP,  CPI,  CP2 
DT 

DTEMP 

DTMIN 

DH 

DHl 

DH2 

DX 

DM 

C0-C3 

DMW 

G 

GTB 

GAMA,  GAMA1.GAMA2 

HF 

HTB 

IBLANK 

IFLAG 

IL.  IR 

IMAX 

JMAX 

MWT 

NALL 

NRT 

NST 

NT 

R 

P,P1,P2 

T 


Equivalent  to  SPECIE  (in  common) 

Concentration 

Throat  area 

Avagados  number 

Areas 

Area  ratio 

Concentration 

Concentration 

Total  mole  fraction 

Specific  heat  table 

Specific  heat  table 

ATime 

A Temperature  for  convergence  check  & 

extrapolation 

A Temperature  for  convergence 
ah  - enthalpy 

ah  - enthalpy  due  to  chemistry 
ah  - enthalpy  due  to  radiation 
Axial  increment  (mm)  input 
A Mach  number 
Temporary  constants 
AMolecular  weight 
Gibbs  Function 
Gibbs  Function  Table 
Ratio  of  specific  heats 
Enthalpy  of  formation 
Enthalpy  Table 
Blank  field 
Write  control 

Stoichiometric  reaction  coefficients 
Number  of  axial  elements  in  the  expansion 

region 

Number  of  axial  elements  in  the  cavity 

region 

Molecular  weight  table 
Number  of  catalytic  species 
Number  of  reactions 
Number  of  species 
Number  of  temps  in  table 
Gas  constant 
Pressure 
Temperahire 

1 
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and  sums  the  values  to  get  the  total  species  production  rates.  These  values 
are  then  returned  to  CHEMS  which  increments  the  species  concentrations 
by  the  necessary  amount  and  sets  up  for  the  next  step. 

A listing  of  the  current  versions  of  these  codes  is  given  in  Appendix  IV. 
The  codes  have  been  run  with  the  input  data  listed  in  Appendix  V;  the  code 
results  have  been  compared  with  a non-lasing  run  with  the  PSI  code.  There 
was  good  agreement  with  the  PSI  code  results  for  this  case;  however,  the 
present  code  should  be  compared  with  a more  thorough  set  of  calculations 
before  it  is  assumed  to  be  totally  valid. 
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SECTION  5 

SUMMARY  AND  RECOMMENDATIONS 


In  the  present  program,  a set  of  general  purpose  subroutines  has  been 
developed  for  use  in  laser  modeling  codes.  These  subroutines  provide 
easy  access  to  collisional  linewidth  data  for  the  HCl  - H2  laser  system  and 
a general  set  of  finite  rate  chemical  kinetics  routines.  Also,  the  linewidth 
and  Einstein  coefficient  data  relevant  to  the  HCl  - H2  laser  system  have 
been  reviewed  and  revised. 

Depending  on  MERADCOM’s  current  interest  and  priorities,  the  pres- 
ent study  suggests  three  possible  actions: 

1)  continue  the  development  of  state-of-the-art  subroutines  with  the 
implementation  of  finite  element  methods  for  solution  of  the 
coupled  gas  dynamics  and  radiative  flow  equations  in  the  laser 
model, 

2)  incorporate  the  present  subroutines  with  presently  available  radi- 
ative and  fluid  flow  routines  to  evaluate  the  HCl  - H2  laser  exper- 
imental results  and  predict  ultimate  laser  performance,  and 

3)  measure  or  analyze  linewidth  and  Einstein  coefficient  data  for 
other  laser  systems  and  format  these  results  for  use  in  the  codes 
developed  here. 

These  programs  could  be  implemented  either  sequentially  or  simultane- 
ously in  the  coming  fiscal  year.  In  any  event,  pursuance  of  these  pro- 
grams would  contribute  to  MERADCOM's  ability  to  accurately  model  and 
evaluate  high-energy  laser  concepts  and  devices. 
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APPENDIX  1 


SELECTED  LINEWIDTH  CALCULATION  RESULTS 

FOR  300  K 
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C MCoL.^.  liD  ;i  ALf“»IOTnS  ?Ur.  hCiO  ?Z' 


^ ? :■ 


1 ■-  r o l'D  d Y .1  C ■ G 


INPUT  CONiTANTS  APU: 


B£(I)  = 1C.1j63!40  B£(F)  = 10.44D4  70  3£(?)  = 10.440470 

BU(bClO)  = 1.123E-13  «U  (HC-0)  = 1.1C8E-18 

Q(HC10)  = 7.60CZ-26  Q (HC-0)  = 7.600E-26 

B(«IS)  = 3.100E-08 


LINE 

HALF-WIDTH 

CNOSS-SECTIOS 

H 

CH-1 

CK**2 

-1 

2.  169063E-D1 

2. 7886413-14 

-2 

2.451 867Z-31 

3. 1524512-14 

-3 

2. 634363E-01 

3. 367736E-14 

-4 

2.579104E-01 

3. 316044E-14 

-5 

2. 323908E-01 

2. 987929E-14 

-6 

1.951782L-01 

2. 509474E-14 

-7 

1.637262E-01 

2. 105C85E-14 

-8 

1 . 367220E-01 

1 . 75  7 883E-14 

-9 

1. 185068E-01 

1. 523684E-14 

- 10 

1 .043747E-01 

1 . 341 982E-14 

-11 

9. 587073E-02 

1.232644E-14 

-12 

8. 892447E-02 

1 . 143334E-14 

-13 

8. 439416E-02 

1.065086E-14 

-14 

8.040339E-02 

1 . 033775E-14 

-15 

7.703322E-02 

9.904431E-15 

- 1 6 

7. 454288E-02 

9.584246E-15 

-17 

7.  191777E-02 

9. 246720S-15 

-18 

7. 04411  9E-02 

9.056969E-15 

-19 

6. 866419E-02 

6. 823397E-15 

-20 

6.736  400E-02 

8.725513E-15 
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c - 


_ 1 


. K 


c t - ; I. ; 

• 0 E E ■ " 1 7 . 

AS  T 

} n^-1  ?■  - :Zr  i 

' J - j AT 

I N ? i : C 7 

s:an:3 

6E(I)  = 

"^.  33  4 .31  E E ( 

‘■)  = 

10. 13o34:  b£(?)  = 10 

. 440470 

.1  J ( ri  J 2 1 ) 

= 1.1 32E- 

1 3 r.  J 

(i;C-0)  = 1.103E-1 

5 

0(.HC2  1)  - 

7.60CE-2 

6 0 (H 

C-0)  - 7.6GGE-26 

B(-,IN)  = 

3 . 100L-03 

LINE 

HALF- »ID7h 

C H 0 

C£-  1 

-1 

9.  91  :^7o3I-02 

1 

-2 

1. 540^90E-01 

1 

- 3 

1 . 399596Z-01 

1 

-4 

1.663214E-01 

c 

-5 

1 . 395C81E-01 

1 

-6 

1. 732100E-01 

2 

-7 

1 . 590666E-01 

2 

-8 

1 . 223545E-01 

1 

-9 

1 .052S46E-01 

1 

- 10 

1.014142E-01 

1 

- 1 1 

9. 1 73584E-02 

1 

- 12 

9. 349430E-02 

1 

-13 

8. 6678o6E-02 

1 

-14 

8. 333015E-02 

1 

-15 

7. 831b99E-02 

1 

- 16 

7. 430738E-02 

9 

- 1 7 

7. 340944E-02 

9 

- 1 6 

7.104^0GE-02 

9 

- 1 9 

6 . 9787091-02 

8 

-20 

6. 829933E-02 

0 

5 -SECT ION 

Cli**  2 

2^b■4^°Z-^U 
9'3D40UE-14 
799509E-1 4 
1 3S452E-  1 4 
793705E-14 
2<;7022E-14 
045176E-14 
573155E-14 
353682E-14 
303318E-14 
1 7943CE-14 
202089Z-14 
114459E-14 
071 4C5E-  1 4 
00694  9E-1  4 
553962S-15 
43851 2E-1 5 
134125E-15 
972769E-15 
7o155CE-15 
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A LC  j 1 A I i.  J .'.A  Li'  - (V  j r H 5 : ^ t i'C  V 2 .-  I l'  F j.  D V iiC  - 0 A I 'F  2 ' . ? i.  - A . K 

I 1>  L.  2 C I ’.’  i r .-.  A 2 S A : 2 : 


£>i,(I)  = 'A.53ybU7  u£(f)  - S,t:3'-I431  dc.(?)  = 10.  '4  40470 

tLMn232)  ^ l.lolE-ld  .“.U(iiC-G)  = 1.1C3i:-13 

a(HC32)  = 7.60CZ-26  Q(riC-O)  ^ 7.60CE-26 

fl(rtlN)  = 3.10Di:-08 


LINE 

HALF-WIDTH 

C.^OES-St-C^I  ON 

CK  - 1 

c a:  « ♦ 2 

-1 

1 . 167527E-01 

1.5C1131E-14 

-2 

1 . 991  6ob£-0 1 

2. 5o0754E-  14 

-3 

1 . 949963Z-01 

2. 5071 36E-14 

-U 

2.099859E-01 

2.  6994  76£-  1 4 

-5 

1 . 780857E-01 

2 . 289710E-14 

-6 

1. 594561E-01 

2.050183E-14 

-7 

1 . 341  285E-01 

1 . 724536S-14 

-8 

1 . 232594E-01 

1. 58479CE-14 

-9 

1 . 285 7 71 £-01 

1.6145892-14 

-10 

1 . 068  U6E-01 

1. 373327E-14 

- 1 1 

1 .030737E-01 

1 . 3252542-14 

-12 

9.45/182E-02 

1.215301E-14 

- 1 3 

3. 929789E-02 

1.143134E-14 

- 14 

8. 422303E-02 

1.032885E-14 

- 1 5 

8 .075905E-02 

1 .033348E-14 

- 16 

7.7  17  162E-02 

9.922232E-15 

-17 

7. 497239E-02 

9. 639465E-1 5 

-18 

7. 210869E-02 

9. 2712702-15 

- 1 9 

7.042955E-02 

9.11 96582-15 

-20 

6. 846426E-02 

8. 60269  5E- 1 5 
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..'MjjLA'rj  ;i  « i HS  ?oi  1'C'43  3y  hc-0  a:  3-.,;,  h 

I !.  ? u r c 0 N o : A r s a r f.  : 


BZ(l)  = 9. 

MU(HC43)  = 
Q(HC43)  = 
B(tUS)  = 


’ 6362  6E(F)  = 9.534647  6E(?}  ^ 10.  4404  70 

1.20SZ-18  MU  (HC-0)  ^ 1.103S-13 

7.faOOL-26  Q(hC-O)  = 7.600E-26 

3. 100E-08 


LI  NE 

H ALF-kIDTH 

csoss-sfc  rioN 

« 

CE-  1 

CM*  *2 

-1 

1 . 494796E-01 

1.921911Z-14 

-2 

1 . 8S373SE-01 

2. 3&341 3E-14 

-3 

2- 2C6487E-01 

2. 3 36  95  3E  - 14 

-4 

2. 043905E-01 

2.62 1920E-14 

-5 

1.401  108E-31 

1.8014S4E-14 

-6 

1 . SSI 008E-01 

1 . 9941 86E-14 

-7 

1 . 55o926E-01 

2. 001  794E-  14 

-8 

1.218  ;2faE-01 

1 . Sob  9S9E-1 4 

-9 

1 . 229S90E-01 

1.S30927E-14 

-10 

1.1413  /9E-01 

1.467S12£'14 

-11 

1 . 037613E-01 

1. 334096E-14 

-12 

9. 304404E-02 

1.260S87E-14 

-13 

9.079379E-02 

1. 167368E-14 

-14 

8.SS671 2E-02 

1.100167E-14 

-IS 

8. 378667E-02 

1.07727SE-14 

-16 

7. 860649E-02 

1 .010672E-14 

-17 

7. 48311 3E-02 

9.6213CSE-1S 

-18 

7.381898E-02 

9.491167E-1S 

-19 

7. 184279E-02 

9.237081E-1S 

-20 

7.043844E-02 

9.CSfaS2'i£-15 
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ISP'oT  CwN 

b-(I)  = 

h J (HC54) 
Q(hC54)  = 

B (;US)  = 


^lLy7!i3 

,-J?  :V.-r-'L  . 'i  ',:C 

- 0 AT  3 0 . 1 L 4 

3 

384  B£(r) 

- 9.23o  36  2 bi.(?)  = 10.4 

^24  70 

1 . 2 33E-  18 

KJ (HC-0)  ^ 1. 108E-1 8 

.630E-26 

Q(iiC-O)  = 7.bOOE-26 

100E-D8 

LINE 

HALF- ki DIH 

CFOSS-5ECriON 

M 

CH-  1 

CX**2 

-1 

1.  43o1d0i--01 

1 . 84654SE- 1 4 

-2 

1 . 659Co4E-01 

2.  133116Z-14 

-3 

2 . 04':#  3 ;6E-01 

2. 634954E-  1 4 

-4 

2. OOuOOLE-OI 

2. 571475E-14 

-5 

1.654608E-01 

2. 127388E-14 

-6 

1 . 835b95E-01 

2. 360474E-14 

-7 

1 . 799435E-01 

2. 313622E-14 

-8 

1 . 474330E-01 

1 . 895593E-14 

-9 

1.020 102E-01 

1.  311561E-14 

- 10 

9.51135&E-02 

1 . 222909E-14 

- 1 1 

1 .027514E-01 

1.  321112E-14 

- 12 

9. 999084E-02 

1 . 28561  8E-14 

- 1 3 

9. 029698E-02 

1. 160980E-14 

- 14 

8. 364051E-02 

1 . 1 3968  2E-14 

- 15 

8. 141C94E-02 

1.046729E-14 

-16 

7. 72041 7E-02 

9.  '426416E-1  5 

-17 

7. 753527E-02 

9.968981E-15 

-18 

7.571 805E-02 

9. 735339E-1  5 

- 19 

7. 219434Z-02 

9. 2822 8 2S -15 

-eO 

6. 722885E-02 

P.6433S3E-15 
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C A L 1 ; 1 f - r.  I ■ .1 S r J I f V.  b b P ~ r I V . r i)  r Y J - J i T i J . ? D . > 

I N P Li  r c ^ I r S A r r. ; 

3E(I)  - b , f'  ? 6 i 7 D bi(r)  = Lr,(r)  = ’'?.4-(jy7C 

MU(nC65)  ^ 1.256E-13  >'.J(riC-0)  = 1.1C8Z-13 

U(HCo5)  = 7.6JCE-26  Q(HC-C)  = 7.6C0E-26 

B(MS)  = 3,1001-33 


LISE 

HALF- WIDTH 

CFOSS-SECriON 

CK-1 

CM**2 

-1 

1. 58o179E-01 

2. 039405E- 14 

- 2 

1 . S38o02E-01 

1 , 973492E-14 

-3 

1 . 921 blOE-01 

2. 470SS3E-14 

-4 

1 . 492031  E-01 

1 . 91 33S7E-14 

-5 

1 . 925095E-01 

2. 47S163E-14 

-6 

1 . 7S9C33E-01 

2.2o16S0E-14 

-7 

1. 85H5^1E-01 

2.  339565E-  14 

-e 

1 . 6C91 1 OE-01 

2,0be339E-14 

-9 

1 . 4237<^8E-01 

1.830537E-14 

- 10 

1 . 2S2192E-01 

1 . 609983E-14 

-11 

1 .074646E-01 

1 . 381  967E-  14 

-12 

1 . 01S50oE-01 

1 . 30S672E-14 

- 13 
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H A 1 E - « 1 0 U n 

U?J35-52Cri'J'i 

W 

C ■!-  1 

C 2 

- 1 

1 .b  1 :(4  32  i,-01 

z.:o2ibi--i- 

-2 

1 . 5 7 1 0 9 1 - : 1 

. 3 2 :■  1 5 b 2 - 1 . 

-3 

1 . 4 4 / 0 ,<  b z - : 1 

1 . 3 5 4 ■'  5 4 b - 1 4 

1 . 5ci1dU3  - )1 

2 . 04 1 94 12-1 4 

-5 

1 .2\,jo-tbi,  — ^1 

1. 5437152-14 

-6 

1 . 6 9d :2  4^- 0 1 

2. 1932162-14 

-7 

1.0311o2£-01 

1.3259C12-14 

-9 

1 . 2 1 u 4 0 1 r.  - 0 1 

1.5b1397i-l4 

-9 

1.492303E-01 

1.93 50492-1 4 

- 10 

1.310365Z-01 

1.6H47332-14 

-11 

1.072319E-01 

1. 37871SE-14 

- 12 

9. 4791712-02 

1.2137712-14 

-13 

7.98375oE- 32 

1.0265022-14 

- 14 

3. 52o62^£-02 

1 .0963232-14 

- 15 

7. 351774b- 32 

9.4524372-15 

- 16 

7. 531 5J6b-02 

9.6935652-13 

- 17 

7 . 376  3 3 1 2-02 

9 .4640032- 1 5 

- 1 9 

7. 1o7450£-02 

9.24  1 16  22-1  5 

- 1 9 

6. 9404362-32 

3 . 9 2 3 56  3 2 - 1 5 

-20 

6. 5145b 12-02 

3. 3760042-15 
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L 


CAL-ULATiJ  HALF- *1370 5 

rOS  i.C6  7 PEi-.Eut-.ijZj  r>Y  HC 

- 1 A : 3 3 3.3  D - ; 

I N ? c : c 0 s s 7 A N : s a? e: 

o£(lj  = 8.D26Ld7  b£(?) 

= 3.  3 350^3  3-  (2)  = K . 1 

d6343 

.Vj(nL87)  = 1.292Z-18 

.1  ;j  ( H L - 1 ) = 1 . 1 3 3 ",  - 1 d 

Q (HC87'  = 7.  bJC£-2  o 

2(h.'-i)  = 7.6::_-23 

3 (MIS)  = 3.10  0r-3d 

L I SI 

H A L : 7 : r. 

c 3 j - j 1 - r : J •' 

a 

2 M - 1 

- 1 

2 . 3c::--?--  3 1 

4.  . 6 D ^ . -1- 

-2 

1 , 7 - i r ^ : 1 

^ . 2 / 3 - : L -1- 

-3 

2 . : . 2 4 3 I - . 1 

2.5  134-;  ,„-1, 

— *♦ 

1 . 7-io  )1 

3 1 3 45'7--1- 

-5 

1 , 2''S71  41-31 

1 . 554  3 ocI-  1 4 

-6 

1.6434551-01' 

2 . 1 1 3 3 5 2 I - 1 - 

-7 

1. 123964--31 

1 . 4451452-14 

- 3 

1.71o7ji2r-)1 

2.237^031-1- 

-9 

1.144  5331- 31 

1 . 4 7 2 '4  ^ - 1 4 

- 10 

1 . 0 1 7 1 3 I - 3 1 

1 . 4?  3o5  5.;-l4 

- 1 1 

1 .0  2d  0711-31 

1 . 3 1 6 7 4 11  - 1-- 

- 12 

9.93ov55I-32 

1 .2o33422-l4 

- 1 3 

3. 3335461-32 

1 . 1 364541-  1 4 

- 14 

3 . 6 1 1 03  5l- 32 

1.1071531-14 

- 15 

7.7203dlE-02 

9.92636oI-1d 

- 16 

3.0134)61-32 

1 . 030  356  I-  1 - 

- 17 

7.7  ^o243I-32 

9.9333:12-13 

- 13 

7. 1454321-32 

9. 137 163E-15 

- 19 

6.39o618E-32 

0. 8674061-15 

-20 

6. 794387E-02 

8. 735735S-15 

w. 

r 


51 


L A 'JL  AT  e:>  half -kid: 

HS 

FOF  fcC93 

P£SrUrt32D  BY 

hC-1  AT  300.0  D20 

INFJ?  co'.'sr  i‘s  a:i*: 

32(11  = 1 . To  6 3^1 1 aZ(?) 

^ 3.C23 

?37  3 3 (p)  = i; 

. 1 j o 3 4 : 

.1  j ( h J 9 3 ) = 1 . 3 J A _ - 

1 3 

3 3 (rX-1) 

= 1.  1 J31 - 1 

d 

^(;iv:9  3)  = 7,o)Cr-2 

6 

0(nC-lJ  = 

7. o?01-2b 

d(^.  iS)  - 3.1D33*j6 

LIME 

n A 1 ■ ■ » 1 7 2 h 

Cr033-3£C71  14 

3 

C:1  - 1 

C.'1=“‘2 

-1 

2.O179D3Z-01 

2. 5945543-14 

-2 

1.54o3o2Z-31 

1. 9339543-14 

-3 

1 . 39J^d03  - 11 

1.7  71351i-l4 

-4 

1 . S7o442E- 31 

2. 4 1 5 1 733-  1 

- 5 

1 . 87^  3743- 01 

2.40  302  13-  1 4 

-6 

1.3477133-31 

1 . 7 3 2 3 : 3 3 - 1 4 

-7 

2 . : 3 D J C 3 - 3 1 

2.5743513-14 

-3 

1.7150453-31 

2. 2025223-14 

-9 

1,  64471  U-31 

2 . 1 1 4 5 6 2 3 - 1 4 

-10 

1. 3o174l:-01 

1 . 7 5 0 3 3 9 - 1 4 

- 1 1 

1 . 275^033-31 

1.639o573-1. 

- 12 

1. 12o135E-31 

1.4479123-14 

- 13 

9 . 577b5o3- 32 

1. 2314333-14 

- 14 

9. 083909E-02 

1 . 1679503-  14 

- 15 

8. 7455213-02 

1 . 1244423-14 

- 16 

8. 0873013-02 

1.0398133-14 

- 1 7 

7. 9541563-02 

1.0226943-14 

- 13 

6 . 3:48603-02 

6.7492473-15 

- 19 

7 . 0 : 1 2 ^33-03 

9.001dG  73-13 

-2C 

7 .0426533-  02 

9. 0549963-  1 5 
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CA  i. ’'JL  A : £0  HALF-«I0TriS 

FOR  HC09  PE?rilcD'D  33  HC-1  AT  309.  9 'J£  3 

INPUT  CJSSTANTS  ARE: 

BE  (I)  = 7.  372670  BE  (F) 

MU(HC09)  = 1 . 3 08E-13 

0(iiCC9)  = 7.600E-26 

= 7.706357  BE  (P)  = 10. 

.>1U(HC-1)  = 1.133E-18 

Q(HC-1)  = 7.600E-26 

136340 

LINE 

H A L ? - « I P I H 

C r 0 3 3 - 3 " 0 7 9 9 N 

« 

C * - * 2 

- 1 

1.978576E-91 

2 . 54  3 92  5E-  14 

-2 

1 . 5835/7E- 91 

2.9359361-14 

-3 

1.421 3 : 1 

1 . 8274  3C1-  1 . 

-4 

1 . 1 31 J 53E- 9 1 

1.4546211-14 

-5 

1 . 57o4o 51-91 

2. 926917^-14 

-6 

2. 275953Z-91 

2 . 92627  31-  1 4 

-7 

2. 1535v0£-91 

2.76883C£-l4 

-8 

1 . 53337oZ-91 

1 . 97794  41-  1 4 

-9 

1 .231 74BZ-91 

1 . 5337  021  -1  4 

- 1C 

1 . 465':^35E-01 

1.384863^-14 

- 1 1 

1 . 1 92  I89i:-91 

1 . 53233  31-  1 4 

-12 

1 . 18595  1E-01 

1 . 52  43  1 9Z - 1 4 

- 1 3 

1 . 102273Z-91 

1 . 4 172321  - 1 4 

- 14 

9.274715Z-02 

1 . 19243  31-  1 4 

- 15 

9.2451451-92 

1.  15S6-!U-1^ 

- 1b 

3. 434c8  3Z-92 

1 . 0 9 : 9 3 2 1 - 1 4 

- 1 7 

7.474724Z-92 

9.607946Z-1 D 

- 18 

7. 0914*5 Z -92 

9.  1 177201-  1 5 

- 19 

7 . 4006  92Z- 92 

9.5153251-15 

-20 

6. 533713E-92 

8. 4070561-  15 
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I W ? J 

3-  (I) 

.1  j ( H : ■ 
5 (HCU 

5 (r. I ,) 


L!\T-3  -:d:h3  f-js  hci:  ? -:-.r'::x a:- d I't  hc-2  ,\:  :'F k 

c ? s 3 : n N r 5 A 3 1 : 


13o343  3?C^)  = K.aO 

47  0 3r(?)  = 9 

. 33  44  3 1 

1 . 1 2 3F  - 1 ^ .1  j (hC-2) 

= 1.167E-1 

a 

7.63':3-23  3(nC-2)  = 
3 . 1 0 3 3 - : 3 

7. oGOE-25 

• T , r 

3 A - W I J : h 

0 3 C 3 i - 0 ' C ^ 4 

Cf-  1 

c r.  = * 2 

- 1 

2 . 264  956F-.01 

2. 912 1343-14 

-2 

2. 5580673-01 

3. 238  = 973-  1 4 

-3 

2. 710969- -01 

3.4855332-14 

-4 

2. 5 3344  2 1-0  1 

3.3230502-14 

-5 

2.249GC51-01 

2.9916242-14 

-6 

1 . 96877G3-01 

2.  4027432-  1 4 

-7 

1 . 5413843-01 

1 . 981 8123-14 

-B 

1 . 327  47  3 r-C1 

1 . 7067792-  1 4 

- ) 

1 . 1532273-01 

1 . 4327452-14 

- 1 C 

1 . ^46837£-01 

1 . 3460202-  1 4 

- 1 1 

9.64312CE-02 

1.24:4=32-14 

- 12 

3 , 966  36  23-02 

1 . 1 52  8372-  1 4 

- 13 

8. 5088733-02 

1 .0940163-14 

- 1 4 

3. 038294E-02 

1.033512E-14 

- 15 

7. 7242853-02 

9. 931337E-15 

-16 

7. U14895E-02 

9.533596E-15 

- 17 

7.2070733-02 

9.2663953-15 

- 18 

7.0230543-02 

9. 0297833-15 

- 19 

6. 907b30£-02 

8. 8813342-15 

-20 

6. 7995C73-02 

8.74236=3-  1 = 
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CAL 


1 S?C 


BL(I) 
So  (HC 
3 {HC: 
3 ( r : ‘ 


H Ao='- *13:33 

rc:-  Hc?i  ?e=:'j:'313  by 

HC-2  a:  39  3..^  313 

A .»  * o A r*  j- 1 

934431  B3(F) 

^ 1 0 . 1 36  3 4 3 BE  (?)  = 

9. 934431 

1. 152F- 19 

-''-■-■(HC-i)  = 1.167;.- 

1 8 

7.60  3E-26 
3.1031-39 

0(HC-2)  = 7.033E-26 

LI  SI 

iiALF-WIOLh 

C r 0 i 3 - 3 "C  " I U'l 

V 

C .A  - 1 

C S • • 2 

- 1 

1.0l4l40Z-31 

1.3039151-14 

— i 

5. 0>i705ot-32 

6.5405051-  1 3 

-J 

1 . 5610401-31 

2.0070451-14 

-4 

1 . 5083732-31 

1.9599741-14 

-5 

2. 3397712-01 

3 .009  325  1 - 1 4 

-6 

1.9551632-01 

2.5139221-14 

-7 

1 . 4717962-91 

1 .99234:1-14 

- 6 

1.101 9042-01 

1. 4167551-14 

-9 

1 . 1615371-31 

1 . 4 934291  - 1 4 

-10 

9.7219941-32 

1 . 2499^111-  1 4 

- 1 1 

3.588.^131-02 

1.1042171-14 

- 12 

8.3515821-02 

1.0737931-14 

- 13 

7. 5353622-02 

9.6884851-15 

-14 

7. 4642602-02 

9. 597063E-  1 5 

- 15 

6.919199E-02 

8.8962622-15 

-16 

6.8431502-02 

8.7984772-  1 5 

- 17 

6 . 490  3392-02 

8. 3455041-15 

-18 

6. 4136512-02 

8 . 2462552-  1 5 

- 19 

6. 2110332-02 

7.9857391-15 

-20 

6 . 143  3 1 32-02 

7.9051011-15 
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1 N ? " r C 0 5 1 A . j \7  7 : 

Q:-(I)  = 9.=-54b47  3^(r)  = 9,934'; 

SJ(H;)2)  = 1.1611-19  «3(.iC-2) 

2(nl3  2)  = 7,o3r;-25  J{HC-Z)  = 

3iri:)  = 3.i:o-?9 


- 1 1 

1 

-3  1 

-4  2 

-5  1 

-0  1 

-7  1 

- 3 1 

-9  1 

-10  9 

-11  3 

-12  7 

-13  7 

-14  7 

-15  6 

-16  6 

-17  6 

-16  b 

-19  6 

-20  6 


PE- 

rOt-SED  Si  HC-2  a: 

39;.'''  D-  ; 

31 

3E  (?)  = 9.9344  3'' 

1 . 16  71-  1 3 

7.60;E-2r 

1 E - • : D ■;  n 

E90S3-3- 1 :i  -!  i 

CH-  1 

C6.4*2 

. '■'2  3 2 74  9-01 

1.315b6 01-14 

. 001 570E-01 

1 . 2877541-14 

. C28388E-01 

1 . 3222351  - 1 4 

. O' 34  70  3 E- 01 

2 . 6 16  03'4:-  1 4 

. 321 316E-01 

1 . 693-’651-14 

. 5577971-01 

2.0C29141-1  4 

. 7079313-01 

2. 1959451-14 

. 34 J676E-01 

1.7276131-14 

. 14/;910E-01 

1.4694901-14 

. 1 17  502E-02 

1 . M221^?-'i  4 

.7264741-02 

1 . 1 22  25  1 1 - 1 4 

.6130451-02 

9.783  3641-  1 5 

. 121 629E-02 

9. 1565311-15 

.0206581-02 

9.026712E-15 

.5127671-02 

8.3736971-15 

.511 2291-02 

8.3717211-15 

. 164388E-02 

7.9257691-15 

. 2C'3  7 161-02 

7.9763361-15 

.2614241-02 

8.0762461-15 

. 0 30  226E-02 

7.7532721-  1 5 
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CA;.  J'.L » : HALr--::U:;^3  HCa3  Pi?:yF3ZD  BY  HC-2  AT  3'';.^  LZ;.  •■ 

I cz'.’s:\;.  rs  a-i: 


23b3o2  3i{r) 

- 3.534647  3B{?)  = 

9.334431 

1 .2332-13 

r J (iiC-2)  = 1 . 1 o7e- 

16 

7 . 6 3 : 1 - 2 6 

3(rtC-2)  = 7.D33F-25 

3 . 1 3 r - 3 

"I 

K A1  r - • 1 1)  r h 

CSOSS-SIZTI  O'i 

V 

C - 1 

C M » » 2 

- 1 

1 . 04q67JB-31 

1.3460025-1 '4 

-2 

1.014261Z-01 

1 . 30407  1 5 - 1 4 

-3 

1 . 5862361-31 

2.0394905-14 

1 . 9417496-31 

2.4965755-14 

-5 

1 . 4 9 3 58  2 - 3 1 

1 . 9267795-  1 4 

-6 

1 . 437386E-31 

1.8433472-14 

-7 

1. 7137905-01 

2. 2034795-14 

_ 3 

1 .2955  365-01 

1.66578':'E-14 

-3 

1 . 1 36420  5-0  1 

1.451 1355-14 

- 10 

1 . 03224  35-  0 1 

1.3914785-14 

- 1 1 

9 . 71  0 109  5-02 

1 . 24845  35-  1 4 

- 12 

8 . 971 8325-02 

1 . 1535465-14 

-13 

8. 3714255-02 

1.0763435-14 

- 1U 

7. 96883 15-02 

1.0245815-14 

-15 

7. 447392E-02 

9.575379E-15 

- 15 

6 . 96563  IE-02 

8. 955957E-15 

- 17 

6. 8371685-02 

8.8550795-15 

- 1 3 

6 . 645C545-02 

8. 5437845-15 

- 1 ? 

6. 4724635-02 

8. 3218825-15 

. ')  A 

^ j 

6. 3016715-02 

8. 1022305-15 

CAic'Ji  :h:. 

iN?jr  ecus: 

B2(I)  = 3. 

a J (riC54)  = 

j (ciCSa)  = 


HAlF-WinrHS  rO 

R HC54  PErrCroi:.  BY 

:3C-2 

A : 3 3 3 ,''  0 . : 

AS  IS  a:-E: 

938384  BS(r)  = 

9.2  3o363  32{?)  - 

9.8344 

3 1 

1.2  332-18  //j 

(:iC-2)  = 1 . 1 672- 

1 3 

7.633S-26  C(H 
3.103^-0- 

C-2)  ^ 7.6  'O: -26 

L I S 2 

li  A 2 2 - W I :>  2 ri 

?0So-3  U”i  ./■< 

2 6.-1 

c a » = 2 

- 1 

1 . 2336192- 31 

1 . 5 8-^.  155-  - 1 4 

_ •-> 

1.4350222-31 

1 . 8452  565  - 1 4 

-3 

1 .966''^0Z-l^ 

2.5278455-14 

-4 

2 . 2 2 3 8 4 0 2 - 3 1 

2.35':>2  705-l4 

-5 

1 . 92864  4.-  -01 

2.4797265-14 

-6 

1 . 8672822-3  1 

2. 40'^8  2 45-14 

-7 

1 . 52952  3.1-01 

1.9655615-14 

_ 0 

o 

1.3611982-31 

1.7521415-14 

- 9 

1 . 3764542-31 

1 . 76  1 755  5-  1 4 

- 10 

9.3895C12-32 

1 . 27  1 528.-  - 1 4 

- 1 1 

1 .021  15o:.-  31 

1 . 3 1 2 1 1 6 L - 1 4 

-12 

9 . 6 22  3co2 -32 

1 . 1 0862  85  - 1 4 

- 13 

8.5832772-32 

1 . 1 0358  35-1  4 

- 14 

8.0748565-02 

1. 0382135-14 

- 15 

7 . 21 J801E-32 

9. 2750 38 £-15 

- 16 

7. 1139525-02 

9.  1 5308'-  E-  1 5 

- 17 

6 . 5681225-02 

8 . 444  36  IS  - 1 ■) 

- 1 8 

6 . 742  364  5-32 

8.668A6^?r-15 

- 19 

6. 246281n-32 

8. 031  ■'625  - 1 5 

-2C 

6. 3678215-02 

8 . 1 87  3 3 2E  - 15 

58 


CAL:jL\''rO  HALF-WIDTHS  FJL  HCS5  PFTTL'PdZD  BY  HC-2  AI  3C3.:  DZ;. 
INPUT  CONSTANTS  APS: 


B£{I)  - 
'‘■j  (HC65) 
; (.iCo5) 
o(TIN)  - 


8.63HI7  0 BS(?)  = 8.93S3R!i  B£(?)  = 9.334i31 

= 1.25cL-18  r.L'(HC-2)  = 1.1671-13 

• 7. 53  r - -26  '3(HC-2)  = 1 , bZ  - 2‘~ 

3 . 1 : I - : 3 


IN; 

HALr-WIDTn 

C 3CSS-35CTI  J.l 

V 

c.r- 1 

C .'« • » 2 

- 1 

1 . 77966  IE-3  1 

2.2881725-14 

_ -1 

i. 

2. 036  4605-3  1 

2.6826355-14 

-3 

2. 2P2202E-01 

2. 9471645-14 

-u 

2.  186b35?-31 

2. 8114975-14 

-5 

1 .3C  65  93S-31 

2. 451 3735-14 

-e 

1.782153E-31 

2. 2913775-14 

-7 

1 . 598037E-31 

2.0546525-14 

- 3 

1 . 346825S-01 

1.7316615-14 

- D 

1 . 1759255-01 

1 . 5 1 1 92«5-14 

- 1 C 

1 . 1254725-01 

1.4470595-14 

- 1 1 

9.321660E-02 

1. 1935195-14 

-12 

9.6370845-02 

1.2455325-14 

- 13 

8.  1035915-02 

1.0419075-14 

-14 

8 . 3068375-02 

1 . 0680405-14 

-15 

8.019423E-02 

1 .0310865-14 

-16 

7. 4922865-02 

9. 6330965-15 

- 17 

7.2427095-02 

9. 3 122095-15 

-18 

6 . 9700965-02 

8.9617005-15 

-19 

6.2335375-32 

8.0146775-15 

-20 

6.5598435-02 

3.4342265-15 

£9 


'1 


CAL,:’Jl?.?r3  d A1  F- iJI  31  .15  F05.  hC^o  PIPIUddZD  Bl  HC-d 


3 ; 3 . .K 


.s  ? V : 


33S5rANT5  .A: 


5£  (I)  = 8.  3 J5523  B5  (F)  = 

«5(HC7d)  ^ 1.2''hT-18  r/J  (.4^.-2)  = 

Q(HC76)  = 7.6)'^£-2b  3(4C-^)  = 

3 (.■!  I N ) = 3 . 1 3 ) - - 3 3 


8.633^7?  BI{?)  = 9.834431 

1 . 1 67E-  1 3 
7 . t : 0 E - 2 o 


IN? 

;i  .A  1 ? - « I 3 I H 

2 r J 5 5 - 5 ~ 2 ” 1 ) N 

X 

C - 1 

cy*: 

- 1 

1.71boo4E-C1 

2.:<^717o2-l-. 

-2 

1 , 723380E- 31 

2. 2 153 i: 2-  1 4 

-3 

2.073  95  12-  3 1 

2.6665522-  1 4 

■"  u 

2.  189«7d2-01 

2.  315o''2"-1  4 

-5 

1.7111442-31 

2 . 2 0 7 0 : - 1 4 

-6 

1 . 25291  32-31 

1.6139112-14 

« 1 

1 . 6 43074  2-31 

2.  1 12557?--'4 

-3 

1 . 59  96342-  3 1 

2.05377"2-l4 

_ 0 

1 . 053  J62  2-  3 1 

1 . 354  34  52  - 1 4 

- 1C 

1.3412442-31 

1 . 7244  3 4Z-  1 4 

- 1 1 

1 . 1 1 2 03.3  2-3  1 

1 . 4 30097  ?-  1 4 

- 1 5 

9. Of  9 3422-32 

1.2817932-1.. 

- 1 3 

3.8304192-32 

1 . 1 322722-  1 4 

- 14 

8.57191  IE-02 

1.1021212-14 

- 1 5 

8.  147907E-02 

1.0476062-14 

- 16 

7. 8172982-32 

1 .0050082-14 

- 1 7 

7.  126  522E-32 

9.  1628202-  1 5 

- 1 3 

6.635427E-32 

8 . 5 3 1 4 3 7 2 - 1 5 

-19 

6 . 7782942 -32 

8 . 7 1 5 3 0 2 ? - 1 :■> 

-20 

6.6732722-32 

8.5762052-15 

60 


c ' - c J : \ ' 

h.Ai?- ;i  t:  Hs 

Fo:  tiC37  ?=■? ruFo:-::;  by 

3 C - 2 A 7 3 0 "'.:  0 3 

I '■  ? j r c j 

s : \ : i A ■ : ; 

31(1)  = 

3.  :2b' 37  3'3(P) 

= e. 33^323  5F(?)  ^ 

9.334431 

' o (.ir  5?) 

= 1.2  3:-£-1S 

X ^ ( h C - 2 ) ==  1.16  71- 

1 3 

3 ( ii  C 3 7 ) r 

7 . 6 2 6 

(nC-^)  ^ 7.600:--26 

BC.-IM)  = 

J.  1 3 

T r « ; 

.i  A 1 ^ - o I 0 . t. 

0 F 3 S 5 - 3 1 2 7 I J 

C''-  1 

CM*  *2 

- 1 

1.b4915j--'1 

2. 1205765-14 

'-2 

1.561C27F-31 

2.9070695-14 

-3 

1 . 9930971-:  1 

2.569:2  - 1 4 

- 4 

1.6637351-01 

2. 1455515-14 

-5 

1 . 390  3701-  :i 

2. 4176651*14 

-6 

1 . 9971841-01 

2.5549935-14 

-7 

1.  3515c:f:-oi 

1. 7376745-14 

-3 

1 . 817683^-01 

2. 3370665-14 

-9 

1.550C12S-01 

1 . 9929:6  5-  1 4 

- 10 

1 . 379393F-01 

1 . 7735355-14 

- 1 1 

1.146436E-01 

1 . 4740735  - 1 4 

- 12 

1.133290S-01 

1.4571105-14 

- 1 3 

7.6344105-02 

9.880118E-15 

- 1 4 

8. 239663E-02 

1 . 059403S-14 

-15 

3. 12631 3E-02 

1 .0443935-  1 4 

- 16 

7 . 712293E-02 

9. 9159545-15 

-17 

7.6093145-32 

9.7335665-15 

- 1 3 

7. 350 165- -02 

9.4503675-15 

- 1 9 

7. 1204075-02 

9.  1549625-  1 5 

-20 

6.6158395-02 

8.5062645-13 

61 


CAiJ’r_\:LD  hal--«id:hs  fc?  HC9j  pi':- jz:'  zy  hc-2  a:  30c.  ■ 3':, 

IS? 3?  coss:a;;:s  a?z: 


B--{I)  = 7.706  357  BE(F)  = 8.C260  37  3Z(P)  = 9.  934431 

SJ(hC99)  = 1 . 304Z-13  «J(HC-  2)  = 1 . 1 672-  1 3 

0 (9C98)  = 7.60  02-26  Q(hC-2)  ^ 7.60  0 2 -26 

6 (r.IS)  = 3 . 1 00  2-0  3 


LIVE 

ri  A.  1 2 - Vi  I J I h 

^ O D 1 / .'1 

•i 

C .'1  - 1 

C.1*  '2 

- 1 

2.0933232-01 

2.6936592-14 

-2 

1 . 749d  97  2-  0 1 

2.2496472-14 

-3 

2.1194:62-01 

2. 7249952-14 

- 4 

1.°166 062-01 

2 . 4 64  32^2-1  4 

-5 

1 . 230o572  - 01 

1.5322332-14 

-6 

1 . 645  3762-0  1 

2. 1 163:42-14 

-7 

2 . 05o  08  92  - 0 1 

2. 6435 962-14 

-8 

1 . 302047  2-01 

2. 3169552-14 

-9 

1 . 4045442-0  1 

1.9058712-14 

-10 

1 . 234  3922-  0 1 

1.5370892-1  t 

- 1 1 

1 .C6b9CO2-01 

1.3717d1;,-14 

-12 

1 . 102  4342-0  1 

1,4174392-14 

- 13 

9. 3966942-02 

1.271 1672-14 

-14 

8. 87o443E-02 

1. 1412762-14 

- 15 

3.6292032-02 

1. 1094872-14 

- 16 

7 . 895C29E-02 

1 . 0 15G922-  1 4 

- 17 

7. 246  05  32-  02 

9. 31651 92-15 

-18 

7.4025392-02 

9 . 5 1 7 0 9 2 - 1 5 

- 19 

7.1033002-02 

9.  1 329672-  1 j 

-20 

6. 855  1242-02 

8.3133792-IS 

caljjia:z^  ■»  1 3:  :i?  f:;-  hcc  ’ ?:::L'r?-ij  f.y  .-ij-j  a:  30'.:  of;. 

I N ? ■;  r c c :j  s : A N 7 3 a ? z : 

dF(l)  ^ 7.37^67}  BZ(F)  = B^{?)  = 9 . B j 4 3 1 

'1J(HC:9)  ^ 1.3:3Z-13  M'J  (:-iC-2)  = 1.167Z-13 

^(HC33)  = 7.ojCF-26  ^ {nC  - 2)  - 7.oCDt-2o 

i3  (F.i.;)  = 3. 1:3:-,-: 3 

LIN"  ;i.AZ"-»I>l  ri  2-jj3-3'2ri'-/I 

.'5  c n;  - 1 c Y » » : 

2. 633o6lZ-l4 
2.0327:21-1-1 

1 . 7 10UQ  5v-  T j 

1 . 7i933n>:-1-, 
3. 22°752Z- 1 4 
2. 97491CZ-14 
2 . 6 5 7 3 M - 1 4 
2.47651  U-14 
2.2011152-14 
2 . 0 : 1 2 1 5 Z - 1 4 
1 . 750  772S- 1 4 
1 . 42505SZ - 1 4 
1.2978352-14 
1 . 132392E-14 
1 .0578022-14 
9.  9121222-15 
1.0592672-14 
1 . 0 143352  -1  4 
9.0431172-15 
8.5414432-15 
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-1 

2.C4o373" -31 

-2 

1 . 6 1 4 35  3Z- 31 

-3 

. 3 3 3 5o7^-31 

- 4 

1.3532232-01 

• c: 

2.  51 1 5a  32-31 

-6 

2.  3137o22-01 

-7 

2.0671742-31 

-3 

1 . 02bi44z-:i 

-9 

1.7119512-3'' 

- 10 

1 . 556  4 76::.-  31 

- 1 1 

1.3616902-31 

-12 

1.1063582-31 

-13 

1.0354112-31 

- 14 

8.811 2352-02 

-15 

8. 227  21  IE-02 

- 16 

7 , 7093012-32 

-17 

8. 2386082-32 

- 18 

7.3895272-02 

-19 

7 . 0 34202-02 

-20 

6 . 6432362-02 

I 


1 


I !.  ? u : : j 


jMD  = 
r. J (hci j) 
^ (rfCiO) 

B (‘■•IN)  = 


1 ; * L,  r * « i I *i  S 

r 3 E,  i 1 J 1 9 t'  3 r I r f b E D b Y 

H C - 3 A : 3 C : 0 3 9 

5. r 3 ^ r r. : 

1 36343  B2 (F) 

= 10.4494  70  qE(P)  = 

9, 534647 

1 . 1 2 3S-  1 8 

S'J{HC-3)  = 1.195E- 

1 b 

7, 600E-26 

Q(HC-3)  = 7.600E-26 

3.  100E-08 

LIKE 

HALF- WIDTH 

EbOSS-SSCTION 

W 

CK-  1 

CK**2 

- 1 

2.3&7676E-01 

2.967063E-14 

-2 

2. 4967983-01 

3. 21C221E-14 

-3 

2. 503091E-01 

3. 2 1831  IE-14 

-4 

2. 295733E-01 

2.951711E-14 

-5 

1 . 963358E-01 

2. 530787E-14 

-6 

1 .679923E-01 

2. 159942E-14 

-7 

1 . 444316E-01 

1 . 857008E-14 

-3 

1. 292578E-01 

1 . 66191 3E-14 

-9 

1.149 102E-01 

1 . 4 7744  IE-14 

- 10 

1 . 033596E-01 

1. 335360E-14 

- 1 1 

9.433 1 1 5E-02 

1.212844E-14 

-12 

8. 839935E-02 

1. 136582E-14 

- 13 

B.349031E-02 

1 . 072397E-14 

- 14 

7. 954383E-02 

1.022724E-14 

- 15 

7. 695502E-02 

9. 894378E-1 5 

- 16 

7. 397C26E-02 

9.510615E-15 

- 1 7 

7. 213527E-02 

9.274686E-15 

- 18 

7.C41 138E-02 

9.053G44E-15 

- 19 

6. 9108C1E-02 

8. 38546CE-15 

-20 

6.825918E-02 

8. 776322E-15 
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L 


1 


I spur  CON 


oPd)  - 

r.  c ( H 0 2 1 ) 

2(HC2  1)  = 


L?- ^ I j: MO 

-j;  nC/il 

'■>  '■  ;< ; i;  n N - !■  ..  y 

ic-3  a:  jm."  ; 

5 A r : 

4 3 1 ( F) 

= 10.  1 363  40  BE  (?)  = 

9.534647 

1 . 1 52E-  1 3 

«U  (HC-3) 

^ 1.195Z- 

1 8 

, bOOE-26 

0(HC-3)  = 

7.600E-26 

100E-03 

LINE 

HALF- WIDTH 

C3OSS-SECTI0S 

H 

CK-1 

C M * * 2 

-1 

1 . 057C43E-01 

1 . 359C83S-14 

-2 

9. 340292E-02 

1.265201E-14 

-3 

2.171 428E-01 

2. 791881S-14 

-4 

1 . 497781E-01 

1 . 925750E-14 

-5 

2. 3C9414E-01 

2. 969295E-14 

-6 

1 . 787698E-01 

2.298507E-14 

-7 

1. U213C3E-01 

1. 827419E-14 

-8 

1. 377743E-31 

1.771U12E-14 

-9 

1.  IIO0I8E-OI 

1 . 423218S-14 

- 10 

1 . 0''6500S-  01 

1 . 294093S-14 

- 1 1 

9. 195573E-02 

1. 182308E-14 

- 12 

8. 467442E-02 

1 . 088639E-14 

- 1 3 

7. 972648E-02 

1 . C25097S-14 

- 1 4 

7. 4 37  988E-02 

9.627573E-15 

-15 

7. 12o57CE-02 

9. 162884E-15 

- 16 

6. 713635S-02 

8.631957E-15 

-17 

6. 255507E-02 

8.042927E-15 

- 1 8 

6. 406  146S-02 

8. 236609E-15 

- 1 9 

6 . 073537E-02 

7.8089593-15 

-20 

6. 158427E-02 

7.9181053-15 
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I'APoV  C0S5: 

d£(I)  - 9. 

rtU{hC32)  = 
3(HC32)  = 

S(KIN)  = 


nALF-oIDrHS  FOF.  liC’2  AF  I ? ^'1 D oY  iiC-3  AT  3C0.^  1- ; A . K 
A ‘i  r 5 A h E : 


5J46U7  BE(F)  ==  9.63  44  3 1 

1 . 1 blE- 18  KU  (HC-3)  = 
7.690E-26  Q(riC-3)  = 

3 . 100E-38 


BE(P)  = 9.5346  47 

1 . 195E-1 8 
7. 600E-26 


LINE 

HALF- WIDTH 

EROSS-SECTION 

H 

1 

u 

CK**2 

-1 

1 . 042287E-01 

1 . 340 105E-1 4 

-2 

5.  257  10  IE-02 

6. 759239S-1 5 

-3 

1 . 495970E-01 

1 . 923421E-14 

-4 

1 . 7C3359E-01 

2. 190C69B-14 

-5 

1 . 88C555E-01 

2.417895E-14 

-6 

1 . 686349E-01 

2. 168193E-1 4 

- 7 

1.54J578E-01 

7. 982060S-74 

- 3 

1. 232982E-01 

1. 535269E-14 

- 9 

1 . 135527E-01 

1 . 4 59  9S8E  - 1 4 

-10 

1 . C47502E-01 

1 . 34681  IE-  1 4 

- 1 1 

8. 84061 7E-02 

1 . 136b95E-14 

- 12 

8. 569378E-02 

1. 101795E-14 

- 13 

7. 722420E-02 

9. 923988E-15 

- 14 

7. 694310E-02 

9.89235GE-15 

- 15 

7 . 063oC9E-02 

9. 081  93  1 S-  1 5 

- 16 

7 .CC3433E-02 

9.004553E-15 

- 1 7 

6.6C2  19BE-02 

8.488630E-15 

-18 

6. 526017E-02 

8. 390729E-15 

- 19 

6.294781E-02 

8.093427E-15 

-20 

6.  226  145E-02 

8.005173S-15 
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1 


.A:£D  nAL?-WIDlHS  FOP,  hC43  PEP'i’UPaSD  BY  BC-3 


30  C 


ISPir  CO“vB:  A.'J  73  ,\P£; 


&F(I)  = 9.23  6362  bF.(F)  = 9,53464?  BZ(P)  = 9.  5346  4? 

3J(HC43)  = 1.2C6F-18  “2(60-3)  - 1.195E-1B 

0(HC43)  = 7.63CF-26  Q(riC-3)  = ?.6CGE-26 

IN)  = 3.i0  3:.-:-t 


L 1 S ■■ 

B L ? - W 1 j 7 K 

OPOSS-SEO'^ION 

n 

C.“-  1 

C.1**2 

- 1 

1 . 212  390E-01 

1. 5594562-14 

-2 

1 . 199b1?E-01 

1.542390E-14 

-3 

1 , 07 1 G99E-0 1 

1 . 377150E-14 

-4 

9.2'^86C5Z-92 

1. 1955542-14 

-5 

1 . 349538E-01 

1. 7351492-14 

-6 

1 . 324773E-01 

1 . 7033072-14 

-7 

1 .022412E-01 

1. 314552E-14 

-3 

1.005  192S-01 

1.292412E-14 

- 9 

9 . 538424E-02 

1.226383E-14 

- 10 

7. 936383E-02 

1.020409E-14 

- 1 1 

3. 471  406E-02 

1 .0891932-14 

- 12 

3 . 170  71  5? -02 

1.050538E-14 

- 1 3 

7. 997353E-02 

1 .0282492-14 

- 1 4 

7 . 199967S-02 

9.2572502-15 

- 15 

7. 241  8752-02 

9.3111382-15 

-16 

6. 638  1752-02 

8.5349342-15 

- 1? 

6. 7297342-02 

8.6526552-15 

- 13 

6. 2994002-02 

8.0993632-15 

- 19 

6.374800E-02 

8. 1963042-15 

-20 

6.  1005292-02 

7.8436642-15 
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JALr'JLATED  h A LF  * k I DT  HS  FOR  IjCSt*  PEFILFBED  oY  HC-j  AT  30'!'. 0 Dii.  K 
INPUT  CONSTANTS  ABE: 


B£(I)  = 8.  93(3384  BS(F)  = 9.236362  BE(P)  = 9.534647 

R0(HC54)  = 1.233E-18  MO  (HC-3)  = 1.195E-18 

U(HC54)  = 7.600E-26  Q(HC-3)  = 7.600E-26 

B (HIN)  = 3 . 100E-08 


LINE 

HALF- WIDTH 

CSOSS-SECriON 

a 

ca- 1 

CK**2 

-1 

2.412975E-01 

3 . 102447E-14 

-2 

2 . 706  469E-0 1 

3. 479802E-14 

-3 

2.857183E-01 

3. 673581E-14 

-4 

2. 763905E-01 

3. 553649E-14 

-5 

2. 497953E-01 

3. 211706E-14 

-6 

2.  111277E-01 

2. 714543E-14 

-7 

1 . 803469E-01 

2. 3 18784E-14 

-8 

1 . 529379E-01 

1.966377E-14 

-9 

1 .349409E-01 

1 . 734983E-14 

- 10 

1.201 172E-01 

1. 544390E-14 

- 1 1 

1 . C85010E-01 

1 . 395036E-14 

-12 

1 .0  1087OE-01 

1 . 2997 19E-  1 4 

- 1 3 

9. 328359E-02 

1 . 199380E-14 

- 14 

8. 871 600E-02 

1. 140679E-14 

- 15 

8.390218E-02 

1 .078760E-14 

- 16 

8.043534E-02 

1. 034186E-14 

- 17 

7. 749659E-02 

9.964011E-15 

- 18 

7. 484323E-02 

9. 622857E- 15 

- 19 

7. 298553E-02 

9. 3840 10E-15 

-20 

7.111 7b4E-02 

9. 143846E-15 
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r 


I 

1 


^ALJliLATLD  HALF-kIOTHS  F3E  HC65  PEInlUPBED  oY  HC-3  AT  300.0  JiG 

i.spur  cossfA:.  rs  afe: 

bE(I)  = 6.63»970  BS(F)  ^ 8,  933384  B£(P)  = 9.  534647 

KU  (HC65)  = 1.256E-13  HU(HC-3)  = 1.195E-18 

2 (hC65)  = 7.6D0E-26  Q (HC-3)  = 7.600E-26 

D (MN)  = 3 . 100Z-03 


LINE 

HALF- WIDTH 

CROSS-SECTION 

H 

CK-  1 

CK**2 

- 1 

1 . 267186E-01 

1 . 629266E-14 

-2 

1 . 463550E-01 

1 . 881738E-14 

-3 

2.015922E-01 

2. 591  94  IE-14 

-4 

2. 366118E-01 

3. 042200E-14 

-5 

1 . 994089E-01 

2.563870E-14 

-6 

1. 964469E-01 

2. 525813E-14 

-7 

1 . 757196E-01 

2.259288E-14 

-8 

1 . 148951E-01 

1 . 477247E-14 

-9 

1 . 445369E-01 

1 . 85836  IE-14 

-10 

1 . 249 183E-01 

1 . 6061  18E-14 

- 1 1 

1 .016675E-01 

1 . 307 174E-14 

- 12 

8. 521789E-02 

1 .C95676E-14 

- 13 

8. 714026E-0.. 

1 . 120393E-14 

-14 

7. 143337E-02 

9. 184443E-15 

- 15 

7.270527E-02 

9. 347978E-15 

-16 

7. 265532E-D2 

9. 341550E-15 

- 17 

6.681770E-02 

8.590991E-15 

-18 

6. 879067E-32 

8. 844660E-15 

- 19 

6.715C24E-02 

8. 633742E-15 

-20 

6. 477964E-02 

8.  328950E-  15 

CALCULATED 

INHUr  CONST 

&£{I)  = 8. 

RU(HC76)  = 
Q(HC76)  = 
a (MM)  = 


HALF-WIDTHS  rOH  hC76  PEFICTbED  UY  liC-3  AT  3CG.G  Dr,G.  K 
ANTS  ARE: 


335323  BE(F) 

= 8.633970  BE(P)  = 

9. 534647 

1.276E-18 

nU  (HC-3)  = 1 . 1 93E- 

18 

7. oOOE-26 

Q(UC-3)  = 7.60DE-26 

3. 100E-08 

LI  NE 

HALF-WIDTH 

CROSS-SECTION 

11 

CK-  1 

CM**2 

-1 

1.815105E-01 

2. 333744E-14 

-2 

2. 120796E-01 

2.  72678  IE-14 

-3 

2. 260234E-01 

2. 906C62E-14 

-4 

2. 225927E-01 

2.861952E-14 

-5 

2.003007E-01 

2. 575337E-14 

-6 

2. 192670E-01 

2. 819192E-14 

-7 

1 . 693 166E-01 

2. 176963E-14 

-8 

1 . 197516E-01 

1 . 539688E-14 

-9 

1.433131E-01 

1 . 842627E-  14 

- 10 

1 .2352C6E-01 

1 . 583 148E-  1 4 

-11 

9. 086698E-02 

1 . 1 68309E-14 

- 12 

1 .0C3016E-01 

1.289613E-14 

-13 

9. 176010E-02 

1. 179792E-14 

- 14 

8. 147246E-02 

1 . 047520E-14 

-15 

8. 163357E-02 

1 .049591E-14 

- 16 

7 . 712048E-02 

9.915652E-15 

-17 

7. 411957E-02 

9. 529819E-15 

- 18 

7. 146454E-02 

9. 188451E-15 

-19 

6. 851387E-02 

8. 809068E-15 

-20 

6.684273E-02 

8. 594210E-15 
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CALJ)LAi-r3  A LF  - • I I)!  HS  FOR  hC37  PKfl'Lr.EFD  bY  3C-3  AT  3C0.0  K 

INPUT  CONSTANTS  An£: 


cE(I)  = 8.02b087  bE(F)  = 8.3388  23  BE(P)  = 9.534647 

flU(HC87)  = 1.292E-18  HU(HC-3)  = 1.195E-18 

^ (HCb7)  = 7.6DCE-26  Q(HC-3)  = 7.6C0E-26 

B (“  IN)  = 3 . 1 00  E- 0 8 


LINE 

HALF- WIDTH 

CROSS-SECTI ON 

H 

CK-1 

CH**2 

-1 

1 . 745334E-01 

2. 244038E-14 

-2 

1 . 745437E-01 

2. 244233E-14 

-3 

2.  1 1 5253E-01 

2. 719655E-14 

-4 

2. 280044E-01 

2. 931533E-14 

-5 

2.0C7617E-01 

2. 581263E-14 

-6 

1. 452291E-01 

1. 867261E-14 

-7 

2.  124242E-01 

2.731212E-14 

-8 

1 . 853359E-01 

2. 382929E-14 

-9 

1 . 583C22E-01 

2.043C62E-14 

- 10 

1 . 390032E-01 

1 . 78721 3E-14 

-11 

1 . 2369G5E-01 

1 . 590332E  -14 

-12 

1.038538E-01 

1 . 335285F  14 

- 13 

9.097314E-02 

1 . 169674E-14 

- 14 

8.953810E-02 

1. 151223E-14 

- 15 

8.005542E-02 

1.029301E-14 

-16 

7. 331824E-02 

9. 426789E-15 

- 17 

7.49a306E-02 

9. 640838E-15 

-18 

7.311431E-02 

9. 4C0629E-15 

- 19 

6.971538E-02 

8. 963550E-15 

-20 

6. 784314E-02 

8.722833E-15 
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CALL,'HLArFD  HALF-WUIhS  rOR  nC?3  PERICFBlD  BY  HC-3  AT  300.0  BEG.  K 
INPUT  CONSTANIS  ARE: 


&£(I)  = 7.  706357  B2(F)  = 8.026087  BE(P)  = 9.534647 

BU(8C98)  = 1 . 304E-18  BJ  (HC-3)  = 1.195E-18 

Q(HC98)  = 7.600E-26  Q (HC-3)  = 7.600E-26 

B (KIN)  = 3.  100E-08 


LISE 

HALF-BIDTH 

CROSS-SECTION 

H 

CM-  1 

CM**2 

-1 

2.  196  474E-01 

2. 824083E-14 

-2 

1 . 544  1C2E-01 

1 . 985306E-14 

-3 

1.961798E-01 

2. 522352E-14 

-4 

1 . 757512E-01 

2. 259695E-14 

-5 

1 . 949 167E-01 

2.506113E-14 

-6 

1 . 898429E-01 

2. 440877E-14 

-7 

1 . 357660E-01 

1. 745591E-14 

-8 

1 . 63381  5E-01 

2. 10D653E-14 

-9 

1 . 463269E-01 

1.881 376E-14 

- 10 

1 . 455439E-01 

1. 871309E-14 

-11 

1 . 236562E-01 

1 .589890E-14 

- 12 

1 .0601  17E-01 

1 . 363030E-14 

-13 

9. 740901E-02 

1 . 252422E-1  4 

- 14 

8.401 1 14E-02 

1 .080  16  IE-14 

-15 

8. 443749E-02 

1.085642E-14 

- 16 

7. 935089E-02 

1 .020242E-14 

-17 

7. 850480E-02 

1.009364E-14 

- 18 

7. 499379E-02 

9. 6422  16E-  15 

-19 

7. 281631E-02 

9. 362252E-15 

-20 

6.716949E-02 

8.636219E-15 
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c aljjla  r 

INPUT  CON 

3 E ( I)  = 

K U (HT09) 
Q(HCC9)  = 
0 (fllN)  = 


D HALF- -.lUTiiS  1 0 E hCC9  F'EEILEjLD  dY  liC-3  AT  J L 0 . G DEG.  K 
TAN  13  AnE: 


372G70  bE(F) 

= 7.706357  bE(P)  = 

9.534647 

1 . 3 C 8 £ - 1 0 

MU  (HC-3)  = 1 . 1 95E- 

1 8 

7. 6aCE-26 

Q(HC-3)  = 7.600E-26 

3. 100E-09 

LINE 

HALF-'nIDIH 

CROSS-SECTION 

K 

CM-  1 

CB**2 

-1 

2. 124305E-01 

2. 731293E-14 

-2 

1 . 859790E-01 

2.391198E-14 

-3 

1 . 445940E-01 

1 . 859095E-14 

-4 

1 .295903E-01 

1.666188E-14 

-5 

1 . 279038E-01 

1. 644504E-14 

-6 

1 . 624 191E-01 

2.088279E-14 

-7 

2. 129009E-01 

2.737341E-14 

-8 

1 . 372969E-01 

1 . 765275E-14 

-9 

1 . 476516E-01 

1.898408E-14 

- 10 

1 . 492765E-01 

1 . 8036 10E-14 

-11 

1. 242890E-01 

1. 598027E-14 

- 12 

1 . 155732E-01 

1. 485965E-14 

- 13 

1 . 032665E-01 

1. 327733E-14 

- 14 

9.441 053E-92 

1 . 2 13972E-14 

-15 

8. 204 180E-02 

1.054841E-14 

- 16 

8. 328277E-02 

1 .070796E-14 

-17 

7. 358629E-02 

9.461250E-15 

- 18 

7.651263E-02 

9.837502E-15 

- 19 

7. 530636E-02 

9.682407E-15 

-20 

6. 965727E-02 

8. 956079E-15 
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APPENDIX  2 


LEAST-SQUARE  FIT  PARAMETERS  FOR  HCl  LINEWIDTHS 


Hcn 

BY 

HC-0 

200 

.0 

8. 

2C750E-02 

1 . 18922EtOC 

- 

1 . 

58161E-02 

-7. 246^7E-C2 

1. 8811''5-P1 

9. 

11270F-C3 

-8.  69417E-C1 

HC21 

BY 

HC-0 

200 

.0 

8. 

36489E-C2 

5.  5:'084S>00 

- 

2. 

7P863E-C  1 

-2.593''  1 B-T  1 

5. C0610E-C  1 

-7. 

4940  3E-02 

-4. 55572E+00 

HC32 

BY 

HC-C 

2C0 

.0 

7. 

61203E-02 

4. 39938E-C1 

- 

149C1S-C2 

1. 32792B-"2 

1. 28617E-C1 

-3. 

1460 9E-C3 

- 1 .95647E-01 

HCU3 

BY 

HC-0 

200 

.n 

7. 

61 389E-C2 

-8.  132‘'4E-01 

3. 

69709S-C2 

1 . ''9826E-C  1 

- 5.  5941''E-02 

1. 

1764  7F-C2 

9.52872E-01 

HC54 

BY 

HC-0 

200 

. 0 

8. 

54757E-C2 

-2.  4459''B  + 00 

6. 

9851 3E-C2 

2.59706E-''1 

-2. 63077S-01 

3. 

C26085-C2 

2. 41 580E*00 

HC65 

BY 

HC-0 

200 

.0 

4. 

4721 9E-02 

-3.9613lEfOC 

3. 

73638S-01 

8. 72833E-02 

-4. 93358E-01 

1 . 

18371E-C1 

3.69823E4-CC 

HC7S 

BY 

HC-0 

200 

.0 

4. 

54566E-02 

-1. 51897EfC0 

2. 

69380B-C1 

-6.  6830  1B-C2 

-2. 60640E-01 

8. 

0412  8E-02 

1 .55741E4-00 

HC07 

BY 

HC-0 

200 

.C 

8. 

90417E-02 

6. 66238E+C0 

- 

3. 

31 357E-C1 

- 3.  12450B-O1 

5. C5881E-01 

-9. 

4754  3F-C  2 

-5. 5T616E+00 

HC98 

BY 

HC-C 

200 

.0 

9. 

77757E-02 

7.  02  942E  + CC 

- 

4. 

C5384E-P  1 

-2.  550:  4E-''  1 

4. 47262B-01 

-1 . 

^24C4E-C 1 

-5. 82C  OOEfCO 

HCC3 

BY 

HC-f' 

200 

.0 

1. 

384C6E-''  1 

9. 8'’815E  + C0 

- 

7. 

93755S-C  1 

- 1.09  1 022-0  1 

8.  3 1 437E-0  1 

-2. 

37619E-01 

-8. 4S386E+00 

Hcn 

BY 

HC-  1 

2C0 

.0 

8. 

21 744E-02 

1 . 3''  356E>0C 

- 

2. 

18672E-02 

-B. 787182-02 

2. 05770E-01 

9. 

39437E-<^3 

-1  ..'ni6PE  + 0C' 

HC21 

BY 

HC-  1 

2C0 

. r 

6. 

0261 5E-C2 

1.93525E^C0 

3. 

5P097E-C  2 

- 2.  26462E-0  ^ 

3.  0 5941  E-C 

3 372  9E-r'2 

- 1 .5685ue>00 

HC32 

BY 

HC-1 

200 

. 0 

6. 

853C'8E-02 

1. 07024E+00 

1 . 

9':'499S-02 

-9.05843E-02 

2.  02  965E-C1 

7. 

93651F-03 

-7. 8S268E-01 

HC43 

BY 

HC-I 

2C0 

. c 

6. 

457C0F-02 

5.  0558-'E-01 

7. 

29623E-02 

-7. 89800B-02 

1. 46191E-C1 

8. 

■'291  2F-C3 

-2»69854E-01 
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HC54 

BY 

HC-1 

200. C 

7. 

510C  3E-02 

-4.  99810E-01 

5, 

45518E-02 

5.991 17E-02 

-2. 565048-02 

1 . 

23f'Cf)1'-02 

6.  7499OE-01 

HC65 

BY 

HC-  1 

2C0.O 

3. 

235368-32 

-2.  3000CE  + 00 

1. 

09197E-01 

1.794268-01 

-2. 863438-01 

5. 

43755E-02 

2. 31264E>D0 

HC76 

BY 

HO-1 

200. c 

5. 

5534  3E-02 

-3.  64636E  + 00 

3. 

31422E-C1 

8.467728-02 

-5.C5111E-C1 

1 . 

170918-31 

3.48743E+00 

HC87 

BY 

HC-  1 

200. n 

3. 

786758-02 

-2. 62918E+30 

4. 

05650E-0  1 

- 1 . 193138-01 

-4.037678-01 

1. 

3C223F-0  1 

2. 57951E+C0 

HC93 

BY 

HC-1 

200.  C 

9. 

784748-02 

7.  51  799E*f^0 

-4. 

15468E-01 

-2.  9D314E-''  1 

5. 6830''8-01 

-1. 

225468-0  1 

-6. 31 216E^OO 

HCC  ? 

BY 

HC-1 

20'!'.  C 

8. 

9693''F-02 

5.04398S*0^ 

-2. 

541  32E-3  1 

-1 .787 18 8-31 

3. 237378-01 

-6. 

95771 8-02 

-a. 15089E+CC 

HC13 

BY 

HC-2 

200. 0 

8. 

383C6F-C2 

1 . 41  365E+-30 

-3. 

34964E-C2 

-7.659738-02 

2. 327048-01 

-2. 

31024E-03 

-1.06534E*00 

HC21 

BY 

HC-2 

200.0 

B. 

7016 1E-02 

3. 16544E*C0 

- 1. 

22231E-0  1 

- 1 . 80879E-01 

4. 382478-01 

-3. 

1 9565E-02 

-2.6643CEtOO 

HC32 

BY 

HC-2 

200.0 

7. 

076078-02 

2.  167roE-t-00 

2. 

43779E-02 

-2. 33835E-01 

3. 16081E-01 

1. 

06769E-02 

-1.76441E+C0 

HC4  } 

BY 

HC-2 

200. 0 

8. 

97669E-02 

3. 19 162E>0C 

-1. 

41 570E-0  1 

- 1.  571  358-01 

3. 440588-01 

-3. 

2861 28-C2 

-2.55799e>00 

HC54 

BY 

HC-2 

200.  F 

8. 

67969F-02 

6. 35936E-01 

-1c 

14624E-02 

1.  17102S--'3 

1. 52666S-01 

-8. 

531278-03 

-3. 68568E-C1 

HC6S 

BY 

HC-2 

200.0 

1. 

0326<"E-01 

-7. 57667E-01 

-5. 

41''33E-02 

2.''5739E-01 

2. 93513S-02 

-7. 

57020’'-03 

7. 3a215E-01 

HC76 

BY 

HC-2 

200.0 

1. 

047888-01 

-2. 51 707e*00 

2. 

53123E-02 

2.  931718-01 

-3. 37550S-01 

4. 

451828-02 

2.54  26iE-^00 

HC87 

BY 

HC-2 

2C0.C 

6. 

14926F-02 

-4.02662E+00 

3. 

51568E-01 

1.^32038-01 

-5.417432-01 

1. 

25  3648-0  1 

3. 86432E^00 

HC98 

BY 

HC-2 

200.8 

7. 

886798-02 

9.  10064E-C  1 

3. 

617578-02 

-9.631428-02 

-1 . 35432E-01 

4. 

1335nf-02 

-4. 32583E-01 

HC09 

BY 

HC-2 

200.0 

7. 

75531E-02 

- 1.27198E+00 

1. 

28975E-C1 

7. 53629E-02 

-3. 15676E-01 

4. 

B2797F-02 

1.45241B*00 

HC10 

BY 

HC-3 

200.0 

8. 

BC781E-02 

9.46121E-01 

-3. 

760078-02 

-1. 914938-02 

1.650668-01 

-8. 

927678-0  3 

-6. 25003E-01 

75 


HC21  BT  HC-3  200.0 

8.60739E-02  1.3063«E>00  -3.52856E-02  -7.57409E-02  2.543025-01 

5.03882E-04  -1.D2503E*C0 

HC32  BY  HC-3  200.0 

8.05356E-02  1.93537E+00  -2.51241E-02  -1.51282E-01  2.82262E-01 

6.43881E-03  -1.54640E+00 

HC43  BY  HC-3  200.0 

7.60850E-02  1.5C503E+00  1.92908E-C2  -1.43752B-01  2.19211E-01 

1.47049E-02  -1.12195E+00 

HC54  BY  HC-3  200.0 

9.21440E-02  1.83172E+00  -8.02218E-02  -8.21555E-02  2.41796E-01 

-7.18594E-03  -1.38712E+nO 

HC65  BY  HC-3  200.0 

7.83078E-02  -1.70299E400  1.464715-01  9.665355-02  -6.157645-02 

3.45812E-02  1.67036E+30 

HC75  BY  HC-3  200.0 

8.81079E-02  -6.57594E-0  1 3.  28420B-C2  1.003445-0''  -3.427775-02 

1.704r5F-02  8.54743E-01 

HC87  BY  HC-3  209.0 

7.32585F-02  -2.25812E+0C  1.69379E-01  1.44325E-C1  -2.906635-01 

5.275365-02  2.27477E+C0 

HC93  BY  HC-3  200.0 

5.63583E-C2  -2.69931E+C0  3.49453E-C1  -2.P5751E-02  -3. 622545-01 

1.05601E-01  2.63630E+00 

HC09  BY  HC-3  200.9 

8.73131E-02  1.48413E>0C  -2.05824E-02  -7.746305-02  -7.75C42E-02 

2.00562E-02  -9.44862E-01 

HC13  BY  HC-0  300.0 

6.33687E-C2  2.39023E-01  2.09231E-02  3.P5279E-03  1.95799E-02 

1.174225-02  -9.01176E-02 

HC21  BY  HC-0  3C0.0 

6.25398E-C2  - 1.74506E*00  1.  137C1E-01  1. ■'78945-91  - 1.433765-C1 

2.951R4E-02  1.50748E+00 

HC32  BY  HC-0  300. 0 

5.59229F-C2  -4.15933E-01  1.C3767E-C1  -4. 113395-03  3.159265-02 

1.506335-92  3.3105^5-01 

HC43  BY  HC-0  3C0.0 

7.  2793  1S-G2  1.47219E^C0  -7.  14972E-C2  - 1 . 1 9 24  9 5- '' 2 1 . 798285-01 

- 3.  35694E-0  2 - 1 . 2 8 1 9^1  E ♦O  9 

HC54  BY  HC-C  300.9 

7.36656E-C2  6.42471E-C2  -1.C1587E-C2  6.079505-92  -5.657005-03 

-2.130645-03  -U.4529CE-03 

HC65  BY  HC-0  3C0.C 

7.  54  380F-C  2 4.  244565-02  -2.454185-0  2 1 . 9 92  3 6 ’ - 9 •'  -4.428145-02 

-1.71201F-C2  6.36592E-02 

HC76  BY  HC-0  3C0.C 

5.83116F-92  1.26012E»00  4.76355E-C2  -1.33362E-01  4.494955-02 

-1.7P371E-C  3 - 9.  42827E-01 

HC87  BY  HC-C  300.0 

6.333255-02  -2.13313E-02  8.32713E-02  -2.9a633E-02  -8.38526E-02 

8.55152E-03  1.88982B-01 
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f 


1 


HC9B  BY  HC-0  3C0.0 


6. 

49040E-02 

- 1. 91044E^0C 

1.62641B-0 1 

1. 35169E-0 1 

-2. 58872E-C1 

2. 

28672E-02 

1. 81 543E+00 

HC09 

BY 

HC-0 

300.0 

9. 

26454E-02 

7.  26205E*r0 

-3. 93217E-C  1 

-2. 76739E-01 

5.05136E-01 

- 1. 

34474E-C1 

-6. 12  168E-*-00 

Hcn 

BY 

HC-  1 

3C0.0 

6. 

35986E-02 

2. 38581E-C1 

1. 98932E-C2 

4. 8101  2S-C3 

1 . 50646E-02 

1. 

32406E-02 

-8.52  98SE-02 

HC21 

BY 

HC-1 

3CO.O 

5. 

93948E-02 

2.  33397Et-G0 

-8. 22282E-C2 

-1. 33961E-"1 

1.6199CE-01 

-1. 

909C8E-02 

-1.99336E*C0 

HC32 

BY 

HC-  1 

3C0.O 

6. 

25222E-02 

3.  26672E-C  1 

3. 36470E-C2 

-1.'>46  13E-C2 

3.  56  950E-02 

8. 

96091E-C  3 

-1  . 67019E-01 

HC4  i 

BY 

HC-1 

300. 0 

5. 

16265E-02 

-9.  86648E-0  1 

1. 0'7278B-0  1 

6.  258  33S-02 

-2. 85647E-02 

2. 

C551 1E-C2 

8.D9068E-21 

HC5U 

BY 

HC-  1 

3C0.  r 

6. 

47129E-C2 

-2.  47202E-C  1 

7. 01425E-03 

9. 16871E-02 

6. 03127B-03 

- 1. 

4C183E-03 

2. 4333SE-01 

HC65 

BY 

HC-  1 

300.0 

6. 

13290E-02 

- 1.  1S244E-C2 

1 . 9C709E-C2 

5. 15 196E-C2 

-5.  3721  3E-02 

4. 

81767E-r'3 

1. 34195E-01 

HC75 

BY 

HC-  1 

300.0 

5. 

89445E-C2 

7. 34663E-C1 

1 . 2254CE-02 

-1. 10 116E-02 

4. 25742E-C2 

-2. 

12974E-02 

-5.53684E-P1 

HC87 

BY 

HC-1 

3C0.0 

7. 

72924E-02 

2. 12833E^00 

-1. 38529E-01 

-1.  18029E-C2 

1. 72174E-01 

-5. 

3712CE-C2 

-1.71 748E+:0 

HC98 

BY 

HC-  1 

300.0 

6. 

54956E-C2 

1 . 5'i046E+00 

-4. 71089E-C2 

1.92523E-C7 

1.36547E-02 

-4. 

77454E-C2 

-1.0'’7‘5SE  + 0C 

HC09 

BY 

HC-1 

300.  C 

5. 

29710E-0  2 

-2. 91 504E+00 

2.332B2E-C1 

1. 79493E-01 

-3.b2192B-01 

U. 

19308P-02 

2. 72  121B1-C0 

HC1  ? 

BY 

HC-2 

300.0 

6. 

20670E-C2 

4, 9a03'iE-C1 

1 . 88381E-02 

-2. 249S4E-r2 

6. 04716S-02 

7. 

860C5E-03 

-3. 1 9268E-C 1 

HC21 

BY 

HC-2 

300.0 

5. 

37690E-C2 

4.95572E-C1 

3. 13628E-02 

-4.S2295E-02 

-1. 31777E-01 

1 . 

98886E-C2 

-2. 71 9^4E-01 

HC32 

BY 

HC-2 

3C0.0 

7. 

'^9730E-C2 

8.  4 1 899B-0  1 

-1.05575E-01 

7. 77603E-07 

-2.69260E-02 

-2. 

32237E-02 

-6.14239E-01 

HC43 

BY 

HC-2 

300.0 

6. 

87D33P-02 

2. D5587B*0C 

-1. 17109E-01 

-3. 25515E-02 

1.20974E-01 

-3. 

96320E-02 

-1.72690B^00 

L 
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1 


HC54 

BJ 

HC-2 

300 

.0 

6. 

07845B-02 

1. 26306B+00 

- 

3. 

98051E-02 

-2.65242E-02 

6. 78234E-02 

-1. 

49747E-02 

-1. 04219B  + 30 

HC65 

BY 

HC-2 

300 

ri 

6. 

11438E-02 

3.02820E-01 

1 . 

93055E-02 

1. 1807BE-02 

3.  3561  IE-02 

-3. 

97976E-04 

-1. 98693E-01 

HC76 

BY 

HC-2 

300 

.0 

6. 

60957E-02 

2. 72001E+00 

- 

1. 

15352E-01 

-8. 99486F-02 

2. 27929E-01 

-5. 

28223E-02 

-2.28  884E-*-3  0 

HC87 

BY 

HC-2 

300 

.c 

7. 

14388E-02 

1. 23904E+00 

- 

7. 

94309E-02 

5.  1 1 896E-02 

5. 92529B-02 

-3. 

9792  1E-02 

-9.70467E-01 

HC93 

BY 

HC-2 

300 

.0 

9. 

149C2E-02 

2. 48447E+00 

- 

2. 

40660S-0  1 

8.  9891  9E-32 

1 . 7949''  S-0’ 

-3. 

29655F-02 

-2. 9"'C66E  + 00 

HC09 

BY 

HC-2 

300 

.0 

5. 

39128E-02 

-2.61  80  3E  *00 

2. 

FF  500E-F  1 

1.98960E-F1 

-4. 04556E-01 

3. 

76474E-02 

2. 53206E+00 

HCID 

BY 

HC-3 

300 

.0 

6. 

53109E-C2 

4. 75639E-01 

- 

6. 

51675E-C4 

2.  199F  1E-C4 

6. 43502E-F2 

-1. 

02016E-03 

-2. 950inE-01 

HC21 

BY 

HC-3 

3C0 

.0 

6. 

n3969E-C2 

1. 0'T719E^CC 

- 

4. 

F6 176E-02 

-1.64803E-02 

2. 588072-02 

-8. 

78938E-03 

-8.72197E-01 

HC32 

BY 

HC-3 

3C0 

.0 

6. 

47326E-02 

2. 03998E^C0 

- 

9. 

78251E-02 

-6. 54136E-02 

2.70146E-02 

-2. 

359C 1E-02 

-1.62762E+00 

HC4  3 

BY 

HC-3 

300 

.0 

4. 

85596F-02 

- 1.63237E+00 

1 . 

59553E-C1 

4. 72764E-02 

-1 . 76932E-01 

3. 

68102F-02 

1. 48675E-t-CC 

HC54 

BY 

HC-3 

300 

. I' 

6. 

42262E-C2 

3.01016E-C1 

3. 

68F99S-C2 

-1. 81293E-03 

2.812902-02 

9. 

2351 4F-03 

-1,34  1 I'^E-OI 

8063 

BY 

HC-3 

3C0 

.F 

6. 

351 12F-02 

1. 4C  087e>00 

- 

5. 

71 250E-C2 

-2.''3039E-02 

6. 96926E-02 

- 1. 

78970E-C2 

-1 .15064E+On 

HC76 

BY 

HC-3 

3C0 

.0 

5. 

86701E-C2 

-6. 31764E-01 

8. 

79095E-C2 

4. 58325E-02 

-6. 57527E-02 

2. 

C0576F-C2 

6. 31 28nE-0  1 

HC87 

BY 

HC-3 

300 

. 0 

8. 

024  1 ''E-C2 

2. 81 48CE+C0 

- 

2. 

C5214F-F 1 

2.  63  76  9E-C2 

2.01 7292-01 

-7. 

56776P-02 

-2. 33427E  + C 0 

HC93 

BY 

HC-3 

3C0 

.0 

6. 

18325F-02 

1 .903C5E*Of 

- 

3. 

51234E-02 

-7. 26800E-F2 

5. 71 332S-02 

-3. 

28486F-02 

-1.44251E+C0 

HC09 

BY 

HC-  3 

300 

.0 

6. 

6D9C4E-02 

-1. 35653E^00 

7. 

48716E-C2 

1. 47456E-01 

- 1 . 39909E-F1 

-8. 

68257F-03 

1 .07864E+00 

Hcn 

BY 

HC-0 

400 

.C 

s. 

4 35‘^0E-02 

1 . 98845E-02 

2. 

36  386F-02 

3.  ''9269E-02 

-2.4 1 485E-C2 

s. 

32499F-03 

7.59366E-02 

78 


HC21 

BY 

HC-C 

400 

.0 

5. 

34164B-02 

1.02781E-C1 

3. 

2P273E-02 

1. 547308-02 

-1. 01 356B-02 

3. 

98248E-03 

-1. 37396E-03 

HC32 

BY 

HC-0 

4C0 

.0 

S. 

9313  1’^-02 

2. 78968E-t-C0 

- 

1 . 

26070E-0  1 

-1. 13495E-C1 

1 . 5795'^E-OI 

-4. 

2^^959^-02 

-2. 33433E+0C 

HC43 

BY 

HC-0 

403 

.0 

6. 

4843  3F-C2 

1 . 258O9E+0C 

- 

7. 

55950E-C2 

1. 457''7E-32 

1. 09C83E-01 

-3. 

39725F-02 

- 1 . 08723E>CO 

HC54 

BY 

HC-0 

4C0 

. C 

5. 

47634F-02 

-3.02808E-01 

5. 

26607E-C2 

4. 83423E-^2 

-1 . 04045E-0  1 

7. 

B3965F-03 

3. 71 8C 1F-01 

HC65 

BY 

HC-0 

4C0 

.0 

5. 

64546E-02 

1.  75863E+C0 

- 

3. 

57133E-02 

-7. 91 399E-C2 

8. 391913-02 

-2. 

56340F-02 

-1.43700E*0C 

HC76 

BY 

HC-C 

400 

.0 

S. 

149Q4E-C2 

-3. 16648E-C1 

1 . 

01 9C4E-0  1 

1 . 66846E-f^2 

-1. 10752E-01 

1. 

26710E-02 

3. 97852E-01 

HC87 

BY 

HC-0 

400 

.0 

6. 

20625E-02 

-3.29181E-02 

1. 

42866F-02 

3.  9445  1S-02 

-5. 93C372-02 

-2. 

G53C7E-02 

1. 1967^E-G1 

HC93 

BY 

HC-'^ 

400 

. 0 

8. 

3070  8S-02 

-3.09482E-02 

- 

7. 

77838E-32 

1 . 84  396E-'i  1 

-6. 12497E-02 

-3. 

99148E-02 

1.538C1E-C1 

HC09 

BY 

HC-C 

4C0 

.0 

5. 

061 95E-02 

- 1.  03005E>00 

1. 

76956E-01 

1. 73 1 09E-O2 

-2.01417E-01 

2. 

54447E-02 

1.03412E*00 

Hcn 

BY 

HC-  1 

400 

.c 

5. 

48654F-C2 

8. 72384E-C3 

2. 

03293E-02 

3. 64266S-02 

-2. 89252E-02 

S. 

62153’='-03 

9.03833E-02 

HC2"' 

BY 

HC-1 

400 

. c 

5. 

51 10OF-C2 

1 . 33645E-02 

2. 

44333E-02 

3.55725E-02 

-2.  77429E-02 

5. 

'^9839E-C3 

8.62126E-C2 

HC32 

BY 

HC-  1 

400 

.0 

4. 

6849CE-C2 

-8.25  362E-01 

1 . 

17579E-C  1 

2. 38906E-02 

- 1. 6944  1E-01 

2. 

8 384  1'^-02 

8.07051E-01 

HC4  3 

BY 

HC-1 

400 

5. 

67451E-C2 

2. 66886E+00 

- 

9. 

1-'410E-C2 

- 1 . 3345  3E-01 

1.  32  1872-01 

-3. 

52754E-02 

-2.22721E^CC 

HC54 

BY 

HC-  1 

400 

.0 

6. 

46308E-''  2 

8. 44  169E-0  1 

- 

4. 

70  34CE-32 

3.63126E-02 

6.  1 4363E-02 

*2. 

51548E-C2 

-7.2f'696E-C  1 

HC65 

BY 

HC-1 

400 

.r' 

5. 

7884 3E-02 

- 1.62178E-C1 

4. 

07230E-02 

5.63C21E-C2 

-1 . 06973B-01 

2. 

6671 9E-03 

2. 60990E-0  1 

HC76 

BY 

HC-1 

4C0 

.C 

5. 

69226E-02 

2. 30297E^CO 

- 

5. 

25581E-02 

-1.07253E-01 

1. 57805E-01 

-4 . 

16968F-02 

-1. 93686E+00 

HC8  7 

BY 

HC-  1 

4C0 

.C 

S. 

86319E-02 

1.00228E+00 

- 

6. 

39406E-02 

5.  96576E-02 

8. 51572E-02 

- 4. 

63765E-02 

-8.5U257E-01 

79 


HC93 

BI 

H2-1 

400 

.0 

6. 

82899E-02 

-5.  41 201E-02 

1. 

09289E-02 

9. 34992E-02 

-6. 852C3S-02 

-2. 

C0848E-C2 

1.65738E-01 

HC09 

BY 

HC-1 

409 

.0 

6. 

2234  IE-02 

- 1.3661  3E4-C0 

1. 

02295E-0  1 

1. 54908E-01 

-2.  16813E-'^1 

9. 

36959E-04 

1. 33437E+00 

Hcn 

BY 

HC-2 

4C9 

.0 

5. 

3P508E-C2 

1. 48114E-01 

2. 

53280E-02 

1. 27799E-92 

-3. 49026E-03 

5. 

40S89E-03 

-3. 83278E-02 

HC21 

BY 

HC-2 

400 

. 0 

5. 

1584  9E-C2 

-1. 3P952E+C0 

1 . 

01 183E-C1 

1.^5951E-'^1 

-2. 75230S-91 

2. 

89427E-02 

1. 35472E+00 

HC32 

BY 

HC-2 

4C  0 

.c 

5. 

99339E-C2 

U.6''127E-01 

- 

2. 

14125E-02 

3. 725423-02 

-1.21 154S-01 

-7. 

9693  IE-03 

-2. 43795E-01 

HC43 

BY 

HC-2 

400 

.0 

5. 

C8536E-02 

- 1 . 29  34  3E  + C0 

1. 

19635E-01 

7. 97C33E-02 

-2. 22692S-01 

2. 

66623E-02 

1 . 235O2E-t-G0 

RCSU 

BY 

HC-2 

400 

. 0 

6. 

14671E-02 

2. 353r4E+CC 

- 

9. 

1 3885E-02 

-9.79732E-02 

9. 32765S-02 

-3. 

13164B-02 

-1.  94C38E4-C0 

HC65 

BY 

HC-2 

400 

.0 

5. 

446C9E-C2 

-2. 4669''E-01 

7. 

15396E-C2 

2. 7841 4E-02 

-5. 87680E-02 

7. 

33329E-03 

2. 43811E-01 

HC76 

BY 

HC-2 

40^ 

.0 

5. 

55762E-02 

8.  3-'56lE-01 

2. 

17843E-C2 

-2. 38037E-O2 

1 .539893-02 

-1. 

71633E-02 

-6.46497E-01 

HC87 

BY 

HC-2 

400 

.0 

6. 

882C7E-02 

2. 5^894E+00 

- 

1. 

21037E-01 

-4.59294E-C2 

2. 20385S-01 

-6. 

95880E-02 

-2. 14592E-4C0 

HC93 

BY 

HC-2 

400 

.0 

7. 

65382E-C2 

1 . 37  376E40C 

- 

1. 

11724E-C1 

9.572973-02 

1.42  464E-''1 

-6. 

76887E-02 

- 1.  19762E  + 00 

HC09 

BY 

HC-2 

400 

.r 

6. 

47761E-02 

-7.06793E-01 

6. 

15685E-02 

1.26159E-01 

-1.3381 92-01 

-1 . 

18931E-02 

7. 37797E-C1 

Hcn 

BY 

HC-  3 

400 

.0 

5. 

46835E-C2 

2.43876E-01 

1. 

C7290S-02 

1. 97 1 1 5E-32 

1. 418712-02 

-1 . 

07372E-03 

-1.21 170E-01 

HC21 

BY 

HC-  3 

400 

.c 

3. 

9431 9E-C2 

4.56 1C2E-0  1 

1. 

03693S-01 

-1. 16139E-:1 

- 1 . 10654S-O1 

2. 

56675E-C2 

-2. 73757E-C1 

HC32 

PY 

HC-  3 

400 

.c 

6. 

10  31  BE-',  2 

-1.32518E+00 

3. 

4^'  389E-02 

1.  4745  33-0'' 

-2. 259823-01 

4. 

63146E-0  3 

1 .03774E*00 

HC43 

BY 

HC-  3 

400 

.0 

7. 

l7678'’-02 

1.  ')9026E400 

- 

1 . 

1 1 785E-C  1 

7. eoc75E-c: 

-3.31 199E-03 

-4. 

6304  1F-02 

-8. 30875E-01 

HC54 

BY 

HC-3 

400 

.0 

3. 

42896E-02 

-1.54227E+C0 

2. 

18714E-0  1 

1. 27938E-02 

-2. 76848S-01 

4. 

71397F-02 

1. 46864E+00 

80 


HC65 

BT 

HC-3 

400.  C 

5. 

76125F-02 

2. 90388E-01 

2. 33915E-02 

3.  13454E-02 

2. 60828E-03 

-8. 

847C7F-03 

-1. 56 180E-01 

HC76 

BY 

HC-  3 

400. C 

6. 

05672E-02 

1. 61  278B-01 

1. 41709E-02 

6. 25860E-02 

-4. 54218E-02 

-1. 

11 174E-02 

-8.73830E-02 

HC87 

BY 

HC-3 

400.  C 

I*. 

1931 1E-02 

4. 80138E-D1 

1. 2‘'486E-C1 

-9. 16294E-02 

-2. 93920E-02 

6. 

968C9E-03 

-3.4D422E-01 

HC99 

BY 

HC-3 

4C0.O 

7. 

7U424F-C  2 

2. 82716E  + 00 

-1.69941E-01 

-3.61212E-02 

2.67C61B-01 

-8. 

r5388E-02 

-2. 43258E4CC 

HC09 

BY 

HC-3 

400. C 

7. 

30037E-02 

8.29926E-01 

-7. 48573E-C2 

1. 1141  IE-01 

5.  95236S-02 

-5. 

63138E-02 

-6. 5871  1E-0  1 

Hcn 

BY 

HC-0 

500.  C 

4. 

77462E-C2 

1. 55645E-02 

2. 54955E-02 

3. 15638E-C2 

-3.04526E-02 

1. 

52178F-C4 

6.  1594  1E-C2 

HC21 

BY 

HC-0 

SCO.  C 

4. 

70922E-02 

5.59751E-C  1 

1. 32655E-C2 

-2. 14774E-02 

1. 84716E-02 

-3. 

43351E-0  3 

-4.63124E-01 

HC32 

BY 

HC-0 

500  . C 

6. 

C5799E-0  2 

2. 36501E^CO 

-1. 42052E-01 

- 4,  67536E-C2 

1. 4C903E-C1 

-4. 

aCS93F-02 

-1.99016E+CC 

HC43 

BY 

HC-C 

500.  C 

4. 

597S3E-02 

7. 33918E-01 

2. 50681E-C2 

-3. 44873E-02 

-2. 62245E-03 

- 1. 

C3217E-02 

-5.  83889E-P1 

HC54 

BY 

HC-0 

500.  C 

5. 

44115E-C2 

4. 29724E-01 

-4. 39796E-03 

3.40090B-02 

4. 00635B-04 

-1. 

40828F-C2 

-3.6044'^E-CI 

HC6S 

BY 

HC-C 

500.0 

5. 

412f''E-C2 

-4.99427E-01 

5. 38519E-C2 

6.75420E-02 

-1. 10353E-O1 

6. 

96327E-03 

4. 91  497E-C  1 

HC76 

BY 

HC-0 

SCO.  C 

4. 

79630E-C2 

9. 44629E-C1 

2. 38627S-02 

-3. 3321 7E-02 

3. 7D577E-02 

-1. 

96779E-f2 

-7. 3C  997E-01 

HC97 

BY 

HC-0 

5C0.0 

5. 

1 873  IE-02 

3.01214E-01 

4. 15168E-C2 

1. 84191E-02 

-1. 49783E-02 

-1. 

80 93 3 F- 02 

-2.066Cf'E-01 

HC99 

BY 

HC-0 

500. C 

4. 

39244E-02 

-5.15063E-01 

1.  26473E-C  1 

2. 41653E-02 

-1. 18731E-01 

2. 

61360E-03 

5.38302E-01 

HC09 

BY 

HC-C 

SCO.  C 

4. 

5517CE-02 

1. 69269E-01 

1.  12C46E-0  1 

- 3.  45586E-02 

-6. 39C27E-02 

-4. 

20359E-03 

- 5.86  524E-02 

HC10 

BY 

HC-1 

500.0 

4. 

82755E-C2 

9.88887E-03 

2. 16061S-02 

3. 93994E-02 

-3. 34400E-02 

-1. 

34450E-04 

7.02531E-02 

HC21 

BY 

HC-1 

500.  C 

4. 

83665E-02 

1.81930E-02 

2. 70933E-02 

3.43617E-02 

-3. 27866E-P2 

-3. 

61147E-04 

6.23419E-02 

81 


HC32  BY  HC-1  500.0 

3.85774B-02  4.31806E-01  6.32056E-02  -5.21512E-02  -6.17334E-02 

5.  66679E-03  - 3 . 0 1 53  8E-0 1 

HC43  BY  HC-1  500. 

6.32320E-02  3. 4642 1 E+00  -2.15224E-01  -6.93958E-C2  2.4DC99H-01 

-7.76987E-02  -2.94907E+00 

HC54  BY  HC-1  500. C 

5.T4939E-02  6.68931E-01  -2.73957E-C3  1.225110-02  2.350220-02 

-1.96473E-02  -5.53139E-01 

HC65  BY  HC-1  500.0 

5.02673E-02  -7.41785E-02  3.18977E-02  6.21762E-02  -4.52902E-02 

-7.85286E-03  8.61451E-C2 

HC76  BY  HC-1  500.0 

4.38696E-02  -1.69187E-01  9.77770E-C2  -1.10386S-C3  -8.72596E-02 

8.25914E-03  2.12707E-01 

HC87  by  HC-1  5C0.C 

5.08926E-02  1.  19024E-t-C0  -7.96675E-C  3 - 2.  33576E-02  5.  831  97E-02 

-3.23471E-02  -9.80202E-01 

HC93  BY  HC-1  500.0 

4.8''914E-C  2 4.  30463E-01  6.57747E-02  - 1.71606E-C2  -1.0491  6E-02 

-1.5340CE-02  -3.01374E-01 

HC09  BY  HC-1  500.0 

5.61397E-02  -2.12714E-01  6.65382E-C2  5.26857E-02  -8.74997E-02 

-1.02222E-02  2.77953E-01 

HCn  BY  HC-2  500.0 

4.  7241  1E-C2  6.53413E-02  2.56518E-'»2  2.4591  8E-02  - 2.  21  396S-02 

1.07333B-03  1.92198E-02 

HC21  BY  HC-2  500.0 

4.77437E-02  3.25049E-C1  -2.28120E-03  3.58716E-02  -4.24C31E-02 

-8.23362E-03  -2.31881E-01 

HC32  BY  HC-2  500.0 

4.  79784F-02  - 4 . 9 0 076 E-0  1 4.00  136E-02  8.  41002E-^2  -1.32071  E-01 

1.23790E-03  5.03971E-01 

HC43  BY  HC-2  500.0 

3.  7551  8E-02  -6.  28680E-01  1.26882E-01  - 1 . 93 1 1 4E-0  ■»  -1.57755E-01 

2.43593E-02  6.24196E-C1 

HC54  BY  HC-2  500.0 

4.47223E-02  1.92863E^00  -1.98441E-32  -1.19255E-ni  7.81696S-02 

-2.32383E-02  -1.614C5E^0C 

HC65  BY  HC-2  500.0 

5.34196F-C2  7.15599E-01  7.35285E-04  1.25134E-02  2.60517E-02 

-2.  1241  9E-'12  -6.02802E-01 

HC76  BY  HC-2  500. C 

4.95398E-02  4.  8652''E-02  4.54207E-02  4.m611E-02  - 4.  32C73E-02 

-8.52O61E-03  -8.34006E-03 

HC87  BY  HC-2  50'^. C 

4.24045E-02  6.4536QE-D1  6.21326S-32  -3.56373E-02  1.59C08E-02 

-1.44942E-02  -5.38141E-01 

HC99  BY  HC-2  50Q.C 

5. 85668E-02  1.64521E*C0  -4.57051E-02  -3.08567B-02  1. 022528-01 

-4.32911E-02  -1.3R180E*00 
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HC09  BY 

HC-2 

500 

.0 

4. 34181E-02 

-6. 89569E-01 

1.48559E-t)1 

2. 67C38E-C2 

-1. 1D247B-P1 

6.  48‘'16F-03 

6. 49585E-01 

HCn  BY 

HC-3 

5C0 

.C 

4. 79058E-C2 

1.57218E-01 

1.62716E-02 

2.  S2479E-'^2 

-6. 82747E-03 

-2. 91680P-03 

-6. 16589E-02 

HC21  BY 

HC-3 

SCO 

.C 

3.7191  1f-'^2 

-4. 33406E-C  1 

1.C26S4E-C1 

-1.71799E-03 

-1.23417S-01 

1 . S8'^92F-02 

4. 4531BE-01 

HC32  BY 

HC-  3 

SCO 

. r 

4. 862  34E-C2 

-1 . 3'^084E  + 00 

6. 67103E-02 

1 . 1SS42E-01 

• 1.  6D207E-01 

6. 51760E-03 

9. 26OO2E-01 

HC43  BY 

HC-3 

SCO 

. r 

5.  ">3310F-?2 

3. 77  307E-C  1 

- 

S.  3341  IE-03 

3. 72167E-D2 

-9.69S20E-02 

-9. 86138E-03 

-2.1 3 109E-F  1 

HC54  BY 

HC-3 

SCO 

. r 

3.  079C  2 

4.  973E-01 

8. 16 179E-C2 

-".80213E-^2 

-9. 86827E-02 

1 . 12S1 3F-02 

-3.  22344E-C  1 

HC6S  BY 

HC-7 

SCO 

. C 

4.  71 96  1F-02 

1.62  623E  + C'' 

- 

9.C7931E-03 

-8.41472E-02 

6. OD481E-02 

-2. 384^9F-02 

-1  . 33  984E  + 0:: 

HC76  BY 

HC-3 

SC'* 

. C 

5. 4934SE-72 

9.  1 2271E-'‘  1 

- 

1.  46  35SE-C  2 

2.  94739E-C  3 

3. 78188E-02 

-2. 46991E-r2 

-7.63242E-C1 

HC8''  BY 

HC-3 

sc  3 

. r 

4.  349B''F-''2 

7. 96144E-01 

6. 4" 181E-C2 

-6. 87873E-02 

4. 61CC0E-03 

-7.  18384F-03 

-6. 19732E-C1 

HC93  BY 

HC-3 

SCO 

.0 

4.  6746  9E-02 

1. 17267E+C0 

3. 61 S29E-02 

-6.2S2S4E-02 

S. 69120E-02 

-2. 5C4C9F-02 

-9. 81  449E-C  1 

HC09  BY 

HC-3 

SCO 

. C 

5. ^23C  3E-02 

5.  46384E-':  1 

5. 27197E-C2 

-4. 17842E-’ 3 

-1. 39539S-02 

-2.  16BS4E-')2 

-3. 8M600E-C  1 

HC10  BY 

H2 

200 

.0 

L-S  FIT  COEPFICIEITS  C (0) 

• • « 

C (HD) 

4. 71013E-02 

1. 241 32E-C1 

1. 12012E-02 

-3. 26742E-C2 

2. 49S85E-02 

-1. 16688E-04 

-4. 32674E-C2 

HC10  BY 

H2 

300 

.0 

L-S  FIT  COEFFICIENTS  C (0) 

• • • 

C (ND) 

3.02429E-02 

1 . 88877E-C1 

- 

4.06S85E-03 

-1 .66012E-02 

2. 40S20E-02 

-3.04313E-03 

-1. 25203E-C1 

HC10  BY 

H2 

40C 

.0 

L-S  FIT  COEFFICIENTS  C(0) 

• • • 

C (ND) 

2. 29643E-02 

1. 38503E-01 

- 

3. 30824E-03 

-9.67S16E-03 

1.74SCC2-02 

-2. 56123E-03 

-1.00183E-C1 

BC10  BY 

H2 

soo 

.0 

L-S  FIT  COEFFICIENTS  C (0) 

• • • 

C (ND) 

1.99680E-02 

7.67147E-02 

5.07801E-04 

-7.65486E-03 

1.00916E-02 

-1.02944B-03 

-5. 70837E-02 
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HC10  BY  AS  2C0.0 

L-S  FIT  COBFPICIEMTS  C(0)...  C (ID) 

1.1435UE-01  5.09231E-01  -3.37538E-01  1.98076E-01  -4.12480E-03 

-2.79543E-02  -2.79080E-01 

HC10  BY  AH  300.0 

L-S  FIT  COEFFICIENTS  C(C),..  C(ND) 

1.31036E-01  -2.10514E-C2  -3.33416E-01  2.28662E-01  -2.29148E-02 

-1.74758E-02  1.13322E-C1 

HC10  BY  AR  400.0 

L-S  FIT  COEFFICIENTS  C (0) . . . C (ND) 

2.25690E-01  -4.41998E-01  -5.14261E-01  3.61613E-01  -4.24690B-02 

-1.69524E-02  4.45504E-C1 

HC10  BY  AF  500. 0 

L-S  Fir  COEFFICIENTS  C (0)  . . . C (ND) 

2.62995E-01  -6.43115E-C1  -5.  82478E-01  4.13800E-01  -5.  44508E- :2 

-1.56157E-02  6.05413E-01 
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APPENDIX  4 

CHEMISTRY  CODE  LISTINGS 
SOBBOOTIME  CBEHS 

C HOFBB'S  CHEHISTBI  CODE,  BODIflED  Blf  TOEB  (3/23/77)  . 

BEAL  HUT 

INTEGEB  BTYP,BCTP,BSPE,B£ACT,ELP,FPP,SPECIE, SPLP,SPFP, AID 
DIHENSION  BE  ACT  (22  5,6)  ,BLP(2  25,6),IPCT(225,6),PRP(225,6), 

♦•HF  (50) 

COHMON/IHEBHO/TTB  ( 2 4 ) , GT E ( 50 , 2 4 ) ,C PT B { 50 , 2U ) , HTB  (50 , 2 U) 

DIHENSION  STE  ( 50,24) 

DATA  IBLANK/'  • / # H 1/82 . 0 5 7/ 

EQUIVALENCE  (STB,GTB) 

DIHENSION  IVT  (5 , 30)  , NVT  (5) 

COHHON/SPECYS/ALPHA  (50)  ,G  (50) 

COHHON/I DENT/ AID  (50)  ,SPLP  (50)  , IS PEC  (50)  , SPRP  (50) 

EQUIVALENCE  (SPECIE, AID) 

COHHON/PATE/BCON  (225,7) ,WDOT (50)  ,RA (225)  ,IPT  (225)  , 

♦ IR(5  0,225)  ,IL(50,225)  ,KB(22  5,3)  ,KL(2  25,3)  , 

♦ aSPE  (225,6)  ,IFLAG  (10) 

COHHON/KONI/CNU  (50) , NST , N ALL , N SPALL , NRT , N T, N E (2)  ,«EXE  (2) 

♦ ,HBT  (50) 

EQUIVALENCE  (IRT,RTYF) 

DIHENSION  PCTP  (2  25)  ,BTYP  (225)  , SPEC  IE  (50) 

C HAIN  PGBAH  FOB  HEL 

C IFLAG(1)=1  PRINT  OUT  THEBHO  DATA 

C IFLAG(2)=1  PRINT  OUT  REACTIONS 

C IPLAG(3)=1  FEINT  OUT  INITIAL  CONCENTRATIONS 

C IFLAG(4)=1  PRINT  OUT  BA  (J)  AND  EV(J) 

C IFLAG(5)=1  PRINT  OUT  BP  (J)  AND  RH  (J) 

READ  (5, 100)  NST,NBT,NT,NALL,  (IFLAG(I)  ,1  = 1,  10) 

100  FOBHAT  (415, 5X, 1011) 

WRITE  (6, 1009)  NST,NBT,NT,NALL,  (I FLAG  (I) ,1=1, 10) 

100  9 FOBH AT  (•  NST=',I5,2X,' NRT=  *,I5,2X,'NT=',I5,2X,'N  ALL=* , 15, 2X, 

1 • NUH= • ,I5,2X, 'IFLAG=' , 1011) 

INPUT  SPECIES  THERMO  DATA 
ONLY  99  SPECIES  CAN  BE  CONSIDERED 
NSPALL=NST>N ALL 
DO  9 NS=1,NSPALL 
DO  ' NR=1,NFT 
IR  (NS,  NR)  =0 
IL  (NS, NB) =C 
9 CONTINUE 

THE  CATALYTIC  SPECIES  WHICH  ARE  DESIGNATED  AS  "ALL"  IN  THE 
PARENTHESIS  (E.G.  HCL  (ALL) ) BUST  BE  THE  LAST  SPECIES  IN  THL 
THERMO  DATA 
DO  10  NS=1,KSPALL 

BEAD  (5, 300)  SPECIE  (NS)  ,SPLP(NS)  ,ISPEC(NS)  ,SPBP(NS)  ,MHT(NS)  ,HF(NS) 

300  FOBMAT(A3,A1,I3,A1,2E10.3) 

DO  11  n=1,NT,2 

READ  (5, 301)  TTB  (IT)  ,CPTB  (NS,  IT)  , STB  (NS,  IT)  , HTB  (NS,IT)  , 

ITTB(IT^I)  ,CPTB  (NS,IT^1)  , STB  (NS , IT» 1 ) ,HTB  (NS , IT* 1 ) 

301  POIHAT  (8P10.  3) 

11  CONTINUE 
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10  COITINOE 

IF (IFL&G  (1) .EQ.O)  GO  TO  13 
WRITE(6,302) 

302  FOHRAT(1H1,‘  THERMO  DATA  FOR  SPECIES  CONSIDERED  IN  THIS  RUN') 

DO  12  MS=1,NSPALL 

HBITE  (6,30  3)  SPECIE(NS),SPLP(NS),ISPEC(NS),SPRP(NS)  ,MHT(NS)  ,HF(NS) 
30  3 FORMAT (II, A 3, A 1,1 3, A 1,2F 10. 3) 

DO  12  IT=1,NT,2 

WRITE(6,304)  TTB(IT),CPTB(NS,IT),STB(NS,IT),HTB(NS,IT), 

1TTB(IT  + 1)  ,CETB(NS,ITt1),STB(NS,IT*1)  ,HTB(NS,IT*1) 

304  FORM  AT ( 11,2  (F10. 0,3F1C . 3) ) 

12  CONTINUE 

13  CONTINUE 

C THE  KL  AND  KF  MATRIX  HAS  THE  COEFFICIENTS  FOR  THE  VARIOUS 

C REACTIONS-NOTE  KL  IS  FOP  SPECIES  ON  THE  LEFT  SIDE  OF  THE  = SIGN 

C AND  THE  KE  IS  FOR  THE  SPECIES  ON  THE  RIGHT  SIDE  OF  THE  = SIGN 

C EXAMPLE  OF  KL  AND  KR 

C KL(J,1)HCL  ♦ KL(J,2)H2  ♦ = KR  (J,  1)  HCL  ♦KR(J,2)H2  ♦ 

C WHERE  J DESIGNATES  THE  REACTION  OF  INTEREST 

C THE  MATRIX  IL  AND  IR  IS  DIMENSIONED  (NUMBER  OF  SPECIES,  NUMBER  OF 

C REACTIONS)  AND  IS  USED  IN  THE  CHEMISTRY  SUBROUTINE  TO  COMPUTE  WDCT 

C THE  ELEMENT  OF  THE  II  AND  IH  MATRIX  WHICH  ARE  NON-ZERO 

C ARE  THE  SPECIES  WHICH  ARE  PRESENT  IN  THAT  PARTICULAR  REACTION. 

C THE  NUMERICAL  VALUE  IS  INDICATIVE  OF  THE  NUMBER  OF  MOLES  PRESENT 

C IN  THAT  PARTICULAR  REACTION. 

C RCTP  IS  THE  PARAMETER  TO  TREAT  REVERSIBLE  OF  NOtl-REVERSIBLE 

C REACTIONS. 

C HTTP  IS  USED  TO  DETERMINE  THE  FORM  OF  THE  RATE  EQUATION 

C RSPE  IS  A MATRIX  WHICH  CONTAINS  THE  NUMBER  OF  THE  SPECIE  RELATIVE 

C TO  THE  ORDER  IN  THE  THERMO  DATA  FOR  EACH  REACTION  WITH  3 REACTANTS 

C AND  3 PRODUCTS 

DO  23  NS=1,NSPALL 
DO  23  IT=1,NT 

HTB  (NS,IT)  = (HTB(NS,I1)*HF(NS))*1000. 
GTB(NS,IT)=HTB(NS,IT)-TTE(IT)*STB(NS,IT) 

23  CONTINUE 

C ♦••••  INPUT  REACTION  RATES  AND  REACTIONS  ♦***♦ 

DO  2600  J=1,NRT 

READ  (5,  14C0)  (KL(J,K)  , PE  ACT  (J,K)  ,FLP(J,K)  ,IRCT  (J,  K)  , 

1 BRP  ( J,K)  ,K=1 ,3)  , (KR  (J,K-3)  , REACT  (J  , K ) , HL  P (J  , K)  ,IRCT  (J  ,K)  , 

2RRP  (J, K)  ,K  = 4,6) 

1400  FORMAT  (6  (II , A3, A 1 ,13  ,A1  , IX)  ,20X) 

8EAD(5,1401)  RCTP(J  ),PTYP(J  ) , (RCON  (J  ,K),K  = 1,7) 

1401  FORMAT  (I1,I2,1X,2F3.1,20X,E10.3,F10.3,E10.3,F10.3,F10.3) 

THE  VECTOR  RSPE  HAS  THE  NUMBER  OF  THE  SPECIE  FROM  THE  ORDER  OF 
THE  THERMO  DATA 

THE  REACTION  RATES  ABE  IN  CM**3/ (PARTICLE*SECONDS) 

DO  2600  K=1,6 

IF(B£iCT(J  ,K)  . ME.  IBLANK)  GO  TO  2400 
BSPB(J  ,K)*-1 


IP(K.LB.3)  KL(J  ,K)=0 
IF(K.GB.4)  KB  (J  ,K-3)*0 

GO  TO  2600 

2400  DO  2500  I*1,NSPALL 

IP(SPECIE(I  ).NE. REACT  (J  , K)  . OB  . IS  PEC  ( I ) . N E . IRCT  (J  , K)  ) GO  TO  250D 
BSPE  (J  ,K)=I 

IP(K.LE.3)  IL  (I  ,J  )=KL(J  ,K) 

IF(K.GE.U)  IB(I  ,J  )=KR(J  ,K-3) 

C GENERATE  IL  AND  IH  IIATRIX 

GO  TO  2600 
2500  CONTINOE 

HRITE  (6 , 1402)  REACT(J,K) 

1402  FOHHAT(1X,'  SPECIE  ',A8,'  NCT  FOUND  IN  DATA  BASE') 

2600  CONTINOE 

IF(IPLAG  (2)  .EQ.O)  GO  TO  2000 
WRITE  (6 ,1403) 

1403  FOBBATCIHI,'  REACTIONS  AND  RATE  CONSTANTS  USED  IN  THIS  CASE  ') 

DO  170  J=1,NBT 

WRITE  (6 , 1404)  J ,(Kl(J  ,K),RSACT(J  ,K)  , RLP  (J  , K)  , IRCT  (J  , K)  , 

♦ RHP  (J,K)  ,K  = 1 ,3)  , 

♦ (KR  (J  ,K-3)  , REACT  (J  , K)  , RL  P ( J,  K)  , I RCT  ( J , K)  , B RP  (J  , K)  , K = 4 , 6)  , 

♦ HCTP  (J  ) ,BTYP  (J  ) 

170  CONTINOE 

1404  FORMAT (1X,IJ,1X,I1,A3,A1,I3,A1,2(*>',I1,A3,A1,I3,A1),'=', 
♦I1,A3,A1,I3,A1,2('^',I1,A3,A1,I3,A1),1X,I2,1X,I2) 

WRITE  (6, 150C) 

1500  FOBMAT(IHI) 

DO  2700  J=1,BBT 
IP  (RTYP  (J)  . EQ.  1) 

♦ WRITE  (6,1501)  J,RCON  (J,1)  ,RCON  (J,3)  ,RCON  (J,4) 

1501  FORMAT  (1X,I3,1X,F3.0,**’,E10.3,'EXP(',F10.3,VT)') 

IP  (RIYP  (J)  . EQ.  2) 

♦ WRITE(6,15C2)  J,RCCN(J,2)  ,RC0N(J,3)  ,RC0N(J,4) 

1 502  FORMAT (1X,I3,1X,F3.0,'*',E10.3,'EXP(',F10.3,'/T)') 

IF  (RTYP  (J)  . EQ.  3)  GO  TC  2700 

IF  (RTYP  (J)  .EQ.4) 

♦ WRITE  (6, 1504)  J,RCON  (J,2)  ,PCON  { J , 3 ) , RCON  ( J , 4 ) 

1 50  4 FORMAT(lX,l3,lX,F3.0,'*',El0.3,'*SXP(',Fl0.3,'*DE23*T*'*1/3)') 

IF  (RTYP  (J)  . EQ.  5) 

♦ WPITE(6,150  5)  J,  RCON  (J,  1)  , RCON  (J,3)  , RCON  (J  , 5 ) , PC  ON  ( J , 6 ) 

1 50  5 FORMAT  (1X,I3,1X,F3.0,'*(',E10.3,'+',E10. 3, ' *EXP (',F10.3,'*DE23*T** 
♦1/3)  •) 

IF  (RTYP  (J)  . EQ.6)  GO  TC  2700 
IP(RTYP  (J)  .EQ.7) 

♦ WRITE  (6, 1507)  J,  RCON  (J,  1)  ,BCON  (J,3)  , RCON  ( J , 4 ) , RCON  (J  , 5)  , RCON  (J  ,6) 

1 507  FORM  AT(1X,I3,1X,F3.0,'*(',E10.3,'*EXP(',F10.3, •♦DE23*T1 3) ♦ ' , ElO. 3, 

♦ ' ♦EXP  ( ' ,F10. 3, ' *T13)  ) • ) 

IP  (RTY P (J)  . EQ.8)  GO  TC  2700 
IF  (RTYP  (J)  . EQ.  9) 

♦ WRITE  (6,  1509)  J,HCOI  ( J , 1 ) ,iCON  (J,2)  , RCON  (J,  3)  ,BCOM  ( J , 5)  ,RCCN  (J  ,6) 
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1509  POHH4T  (1X,I3,1X,  F3.0,'*  • ,F3.0,'*(',E10.3,'<-',B10.3,'*EXP(',P10.J,' 
♦•DB23*T**-1/3) ) • ) 

IF(RITP(J)  .EQ.  10) 

♦ HR  IT  B (6, 1510)  J,RCON  (0,1)  ,BCON  (0,2)  ,RCON  (0, 3)  ,RCON  (0 , 4 ) , RCON  (0 , 5)  , 

♦ PCON  (0,6)  ,BCON  (0,7) 

1510  FORMAT  (1X,I3, IX, F3.0  F3.0, •♦(',E10.3,' 5* E XP  ( ' , F 1 0 . 3 , ' ♦ D E 

♦)>',E10.3,'*T**.5*EXP(-(',F10.3, •♦DE23*' ,F10.3,')*T»*-1/3) ') 

IF (RTY P (0) . Eg. 1 1 ) 

♦ WRITE  (6,1511)  0,RC0N  (0,1)  ,RCON  (0,2)  ,RCON  (0,3)  ,PCOS  (0,4)  ,RCON  (0,5)  , 

♦ PCON  (0,6)  ,RCCN  (0,7) 

1511  FORMAT  (1X,I3,1X,F3.0,'*',F3.0,'*(',E10.3,'*T**-.5*EXP(',F10.3,'*DE 

♦ ♦•,E10.3,'*T**.5*EXP(',F10.3,'  *DE2  3^  • , F 1 0 . 3 , ’ ) 1/ 3)  • ) 

2700  CONTINOE 

2800  CONTINUE 

READ (5, 59  9)  (HE (I ) , HE XE  (I)  , 1= 1 , 2) 

599  FORMAT (4F1C. 8) 

READ  (5,600)  (CNO  ( NS)  , NS= 1 , NST) 

600  FORMAT  (7E1C. 3) 

IP(IFLAG(3)  .EQ.O)GO  TC  18 
WRITE  (6,601) 

601  F0RMAT(1H1,'  INITIAL  CONCENTRATIONS  OF  SPECIES  IN  MOLE  ', 

1 •FRACTIONS  AT  NOZZLE') 

WRITE  (6,602)  (SPECIE(NS),SPLP(NS)  , ISPEC  ( NS ) , SP  RP  ( NS ) , CN  U (K  S)  , N S=  1 , N 31) 

602  FORMAT (3  (2X,A3,A1,I3,A1,  • = ',  1PF11.4) ) 

18  CONTINOE 

C CONVERT  FROM  MOLE  FRACTION  TO  MOLES/CC 

CNTI=0. C 
HTVB=0.0 

C HTVR  IS  IN  UNITS  GRAMS  GAS  /HOLES  GAS 

DO  19  NS=1,NST 
HTVR  = HTVR ♦CNO  (NS) *MWT (NS) 

MUST  HAVE  NON-ZERO  CONCENTRATION  POP  RATE  OF  PRODUCTION  OF  SPECIES 
EQUATIONS  TC  HORN  PROPERLY 
CNTI=CNTI^CNO (NS) 

19  CONTINUE 
HRITE(6,603)  CNTI 

603  FORHAT(1H0,'  THE  TOTAL  MOLE  FRACTION  =',F10.8) 

CALL  CHEHVT  (I  VI, NVT) 

RETURN 

ENTRY  CHEM  (T,P,DT) 

HRITE(6,611)  T,P,DT 

611  F0RHAT(1H1,'  TEMPERATURE  = • , 1 PE  1 1 . 4 , 2X , * PR ES SURE  ( ATM)  =', 

♦ 1PE1  1.4, 2X, 'DELTA  TIME  =',1PE11.4) 

C RHO  IS  IN  UNITS  GRAMS  GAS/CH**3 

BHO  = P*HTVR/  (T*R1) 

HRITE(6,607)  RHO , HT V E, SOM DT 

607  FORM  AT  (2X, 'RHO  = • , 1 P E 1 1 . 4 , 5X , • MOLECU L AR  HEIGHT  OF  GAS  =', 

♦ 1PE1 1. 4,5X, 'SUM  DELTA  TIME= • , 1 PE  1 1 . 4 ) 

DO  20  NS1=1,NST 

IF  (CNO  (NS  1)  .LB. 0.0)  CIO  (IS1) =1.0E-50 
ALPHA  (IS  1) -CIO  (IS  1) /HTVR 
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CM  MOM  IS  CCMT8RT  PBCH  HOLE  PB&CTIOM  TO  HOLES/GH  GAS 

20  ALPHA  (MSI) =ALPHA  (MSI) *BHC 

C CM  AND  ALPHA  ARE  NOW  IN  UNITS  OF  HOLES  I/CH**3 

CALL  CHHSDH  (HT,NVT,  ALPHA  (1)  ) 

NBITE (6,608) 

608  POBHAT(1HO,‘  INITIAL  CCNCENTHATIONS  OF  SPECIES  IN  MOLES/CH ♦ ♦ 3 • ) 
HRITE(6,60  2)  (SPECIE  (NS)  ,SPLP  (NS)  , I SPEC  (NS)  ,SPHP  (NS)  , 

♦ALPHA(NS) ,NS=1,NSPALI) 

CALL  CHHSTY  (T, ALPHA  (1) ) 

NPITE  (6,610) 

610  FORMAT  (1H1) 

NRITE  (6,609) 

609  PORHAT  (2X, 'SPECIE* ,8X, 'C NU ' , 2 1 X , • NDOT ' / 

115X,  • (MOLE',19X, • (MOLES/*/ 

213X, • PH ACTIO NS) * , 1 UX ,* HOLES- SEC)  *) 

DO  21  NS2=1,NST 

NDOT  IS  IN  UNITS  HOLES  1/ (CH**3-SECOND) 

HDOT  IS  NON  CONVERTED  TO  UNITS  MOLES  1/ (HOLES-SEC) 

BOOT  (NS 2)  =HDOT  (NS2)  *MTVR/HHO 
CNU (NS2) =CNU (NS2) +HDCT  (NS2) *DT 
IP(CNO  (NS2) .LT.O.O)  CNU  (NS2) =1 .0E-50 

MBITE (6,605)  SPECIE (NS2)  ,SPLP (NS2)  , I SPEC  (NS2) ,SPRP (NS  2)  , 

♦ CNU  (NS2)  ,HDCT  (NS2) 

605  FORMAT (IX, A3, A1 ,I3,A1,*=*,1PE15.8,10X,1PE15.8) 

21  CONTINUE 
IFITE (6,612) 

RETURN 

612  FORMAT  (1H0,*  CONCENTRATIONS  OF  SPECIES  IN  MOLE  FRACTIONS  AT  DELTA 
NRITE  (6,606)  (SPECIE  (NS)  ,SPLP (NS) , IS  PEC (NS) ,SPPP (NS)  , 

♦CNU(NS) ,NS=1,NSPALL) 

606  FORMAT (3 (1X,A3,A1,I3,A1,*=*,1PE11.4)) 

FND 

SUBROUTINE  CHEMVI  (I VT , NVT) 

INTEGER  SPECIE, SPLP,SPRP, AID 

COMMON/IDENT/AID  (50)  ,SPLP  (50)  , IS PEC  (50 ) , SP B P (50 ) 

COHHON/KON1/CNU (50),NST,NALL,NSPALL,NRT,NT,HE(2),NEXE(2) 

♦ ,HNT  (50) 

THIS  SUBROUTINE  CALCULATES  THE  ARRAYS  USED  TO  COMPUTE 
THE  SUBS  OF  CONCENTRATIONS  IN  THE  VT  DEACTIVATION  PROCESS 
H IS  THE  SPECIE  WHICH  IS  TO  BE  ADDED  TO  THE  "ALL"  COMPOUNDS 
L IS  THE  "ALL"  COMPOUNDS  POSITION  IN  ARRAY  I VT 

NVT  IS  THE  NUMBER  OF  SPECIES  THAT  WILL  CONTRIBUTE  TO  THE  "ALL" 
COMPOUND 

DIMENSION  IVT  (5, 30) , NVT  (5)  , lAI  (8) 

DIMENSION  SPECIE  (1) 

EQUIVALENCE  (SPECIE, AID) 

DATA  ILPAEN/1H  (/,IALL/999/, L/0/ 

NOTE  ALL  SPECIES  WITH  ALL  IN  THE  FORMULA  MUST  BE  AT  THE  END 
OP  THE  THERMODYNAMIC  DATA 
NSP1=NST*  1 
DO  6 I*NSP1 ,NSPALL 
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DECODE  (8, 100, SPECIE  (1)  ) (HI  (K)  ,K=1,8) 

101  POBH&T  (lX,He,8  (IX, k1) ) 

100  F0BRAT(8A1) 

DO  1 11=1,8 
J = I1 

IF(SPLP (I) .EQ.ILPAPN)  GC  10  2 

1 COETINOE 

HRITE(10,1000)  SPECIE(I) 

1000  FORMAT (IX, A3, A1, 13, A1, 'THE  SPECIE  DOES  NOT  HAVE  PROPER  LABEL') 
GO  TO  6 

2 IF(ISPECd)  .EQ.IALL)  GC  TO  3 

WRITE  (10,  1000)  SPECIE(I)  ,SPLP(I)  ,ISPEC(I)  ,SPRP(I) 

THE  LAST  SPECIE  WITH  ALL  IN  ITS  FORMULA  HAS  BEEN  FOUND 
GO  TO  6 

3 COHTINOE 
! = !♦  1 

IF(J.EQ.3)  ENCODE  (2,205, ICHAR)  (I  AI  (K ) , K=  1 , 2) 

205  FORMAT  (2A1) 

200  FORMAT (A2) 

IF(J.EQ.4)  ENCODE  (3, 305, ICHAR)  (I  AI  (K)  , K=  1 , 3) 

305  FORMAT  (3A1) 

300  FORMAT  (A3) 

N=0 

DO  5 M=1,NSI 

IF(J.E0.3)  CECODE(2, 200, SPECIE  (M))  ICHAP1 
IP(J.EQ.U)  DECODE  (3, 300, SPECIE (M)  ) ICHAPI 
IF  (SPECIE  (M)  . EQ.  SPECIE  (I)  ) GO  TO  U 
GO  TO  5 

4 N=l4l 

IVT (L,N) =M 

5 COMTINOE 
Nvr (L) =N 

6 COHTIHOE 

7 RETURN 
END 

SUBROUTINE  CHMS UH  (I VT , N VT ,CN) 

DIMENSION  CN (1) 

COHHON/KON1/CNU (50)  , BST , NA LL , N SP ALL , NRT , NT, HE ( 2)  ,HEXE  (2) 

♦ ,MHT  (50) 

DIMENSION  IVT  (5, 30)  , NVT  (5) 

THIS  SUBROUTINE  COMPOTES  CONCENTRATIONS  FOR  THE  COMPOSITE 

COMPOUNDS 

III=NST* (LL-1) 

DO  2 I=1,NA11 
L=NST*I 
CM(L) =0.0 
N = NVT  (I) 

DO  1 K=1,i 

1 CI(L)  =CN  (1) +01  (IVT  (I,K)  ) 

2 CONTINUE 
RETURN 

END 
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SOBPOOTIIE  CHHSTT  (T,CH) 

«*«*•  SOBPBOGHAH  1 - CHEHISTBr 

*♦*  DNITS  = KILCCALORIES, MOLES, CUBIC  C EN TI METERS  , SECON DS , DEG  RE £S  KELVIN 
***  INPUT  = CN-HIXED  SPECIES  CONCENTRATION  (MOLES/CC) 

IL  AND  IR  APE  THE  STCICHICMETPIC  COEFFICIENTS  FOR  THE  REACTIONS 
IL-LEFT  SIDE  SPECIES/REACTION  CORRELATION 
IR-BIGHT  SIDE  SPECIES/REACTION  CORRELATION 
TP-TRANSLATICNAL/ROTATIONAL  TEMPERATURE 
6 TEMPERATURE  EXPANSION  PARAMETERS  XLGR  PER  ALOG  OF  REACTION  NR  PATE 
COEFFICIENT  RA(NP),  6 TEMPERATURE  EXPANSION  PARAMETERS  XLGK 
PER  SPECIES  NS  EQUILIBBIUM  COEFFICIENT  ALOG  XLGRP(NS), 

6 TEMPERATURE  EXPANSION  PARAMETERS  HFX  PER  SPECIES 
NS  HEAT  Cf  FORHATICN  HF  (NS)  (KCAL/MOLE) 

NST-TOTAL  SPECIES  NUMBER  , NRT-TOTAL  REACTION  NUMBER  *♦* 
OUTPUT  = TEC-CHEMICAL  PRODUCTION  THERMAL  ENERGY  (KCAL/CC/SEC) 

RCN-CHEMICAL  RATE  OF  FORMATION  OF  SPECIES  (MOLES/CC/SEC) 
RCNT-TCTAl  CHEMICAL  FCBMATION  RATE 
PRESCRIPTION  FOR  REACTION  NR= 

PDT  (CN  (NS)  **IL  (NS,  NB)  ) - (RA/RV)  - PDT  (CN  ( NS)  * *IR  ( NS , N P. ) ) 

WHERE  POT-PRODUCT  OVER  ALL  SPECIES  NS  ***♦♦ 

INTEGER  HTYP,RSPE 
REAL  HHT 

COHMON/BATE/BCON  (225,7) ,WDOT(59) ,RA (225) ,IPT (225) , 

♦ IR (50,22  5)  ,IL  (50,225) ,KR  (225,3)  ,KL  (225,3)  , 

♦ BSPE  (225,6) ,IFLAG (10) 

EQUIVALENCE  (IRT,RTYP) 

DIMENSION  BTYP  (225)  ,FV  (225) ,PP  (225) ,RH (225) 

COHHON/SPECYS/ALPHA (50) ,G (50) 

COHMON/KON1/CNU (50) , NST, NALL, N SPALL, NBT, NT, WE (2)  ,WEXE (2) 

♦ ,HWT  (50) 

DATA  AVAGAD,R/6.0228E23, 1.98717/ 

DIMENSION  CN (1) 

BTP=R*T 

TPB=T-298. 

CALL  HATH(T) 

IF(IFLAG  («)  .EQ.O) GO  TC  1 
WHITE(6,100) 

100  FORHAT(1X,'  J • , 1 0X  , » R A ( J)  • , 1 5 X , ' B V ( J)  • ) 

1 CONTINUE 

DO  15  J=1,NET 

YKPT=0.0 

NUL=0 

NUB=0 

DO  5 K=1,3 
KL1=  RSPE(J,K) 

IF(KLI.EQ.-I)  GO  TO  5 
NUL^NUL^KL  (J,K) 

YKPT  = YKPT-KL (J,K) *G  (KL1) 

5 CONTINUE 
DO  10  K=1,3 
KB1  = RSPE  (J,K^3) 
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IF(KB1 .BQ.-1)  GO  TO 
N0R=N0R^KP  (J,K) 

TKPT  = TKPT*KR  (J,K) *G  (KB1) 

10  COBTIRDE 

IMBU>IOL-IOB 

C THE  REUCTION  BATES  ARE  IN  Cfl** 3/ (PARTICLE* SECONDS) 

BA(J  )=RA(J  )*AVAGAD 

C FATES  ARE  NCW  CONVERTED  TO  CH**3/ (BOLE*SECOND) 

RV(J  )=  RA(J  ) /EXP  (-YKPT/RTF) 

IP(RCTP  (J)  .EQ.  1)  PV(J)=0.0 

NOTE  IP  BCTF  =1  THEN  BEACIION  IS  NOT  REVERSIBLE 
IP(INEU.NE.O)  RV(J  )=RV(J  )*RTP**INEU 
IP(IPLAG  (4) .EQ.O)GO  TC  15 
HRITE  (6, 101)  J.RA  (J)  ,BV  (J) 

101  FORHAT  (1X,I3,5X,E15.e,5X,E15.S) 

15  CONTINUE 
ISET=1 
IPRINT=1 

IF(IFLAG  (5)  .EQ.O) GO  TC  16 
WRITE  (6 ,200) 

200  FORnAT(1X,'  O' , 1 0X , • RP  ( J)  • , 1 5X  , • RM  ( J)  • ) 

16  CONTINUE 
DO  108  1=1, NST 
NDOT  (I) =0.0 
DO  107  J=1,NBT 
IF(ISET.EC.C)  GO  TO  25 
RCN1=1.0 
RCN2=1.0 

DO  20  NS1=1,NSPALL 
FCN1=HCN1*CN (NS1)  **Il (NS1,J) 

RCH2=RCN2*CN  (NS1)**IF(NS1,J) 

20  CONTINUE 
BP  (J)  =HCN  1 
BH  (J)  =RCN2 

IP  (IFLAG  (5)  . EQ.O) GO  TC  21 
HPITE(6,2C1)  J,RP  (J)  ,RH  (J) 

201  FOBHAT  (1X,I3,5X,E15.8,5X,E15.0) 

21  CONTINUE 
25  CONTINUE 

BKCN  1 = RP  (J) *RA  (J) 

BKCN2=RM  (J)  *BV  (J) 

HDOT  (I)  =HDCT  (I)  ♦ (IR  (I,  J)  -IL  (I,  J)  ) * (BKCN  1-RKCN2) 

107  CONTINUE 
ISET=0 

100  CONTINUE 

WRITE  (6,1016) 

1016  FOBHATC  FINISHED  CHRSTY') 

RETURN 
END 

SUBROUTINE  HATH  (TP) 

INTEGER  BTTP,BSPE,BEACT,BLP,BRP,SPECIE,SPLP,SPBP 

• « 
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COHHOI/THEBBO/TTB (24) ,GTB  (50,24) ,CPTB  (50,24) ,HTB  (50,24) 
COHBOB/BATE/BCOH (225,7) ,BD0T (50,10) , PA  (225)  ,IBT  (225)  , 

♦ IP (50,225) ,IL  (50,22  5) ,KP  (225,3) ,KL  (225,3) , 

♦ RSPE  (225,6) ,IFLAG  (10) 

EQDIVALEHCE  (IBT,FTIP) 

DIBENSION  PTyP(1) 

COMBOB/SPECIS/ALPHA  (50)  ,G  (50) 

CCMBON/KONI/CNO  (50)  , MST, NALL, N SPALL, NET, NT, 4E  (2)  ,4EXE  (2) 

DATA  JTKEy/1/,R/1. 98717/ 

P23=2./3. 

JP=JTKEy*1 
10  SDT=TP-TTB  (JTKET) 

HDT*TTB  (JF)  -TP 
IP  (SDT*HDT)  15,52,50 
15  If(SDT)  30,40,20 
C POSITIVE  STEP 

20  JTKEy=JP 
Jt=JP* 1 

IP(JTKBy-NT)  10,40,40 
30  JP=JTKEY 

JTKEy=JTKET-1 
C NEGATIVE  STEP 

IP(JTKEY.GE.O)  GO  TO  10 
40  HRITB(6,100)  TP,JTKEY 
C OOT  OF  RANGE 

100  FOBBAT(1H  ,•  TEMPEPATDBE  GOT  OF  PANGE-CHEB  ',Fm.5,I5) 

STOP 

50  SDT»SDT/ (SCT*HDT) 

GO  TO  54 

52  IF(HDT.EQ.0.)SDT=1. 

54  HDT=1.0-SD1 

DO  60  I=1,NSPALL 

G(I  )=GTB(I  ,JTKEY)  ♦HDT^GTB  (I  ,JP)  *SDT 
60  CONTINUE 

HTTP  SPECIFIES  THE  FUNCTION  FORM  OF  THE  REACTION  RATS.  READ  IN 
HHEN  RATES  APE  READ. 

THE  REACTION  RATES  ARE  IN  Cfl**3/ (PARTI  CLE*3ECONDS) 

T13=TP**  (-1./3.) 

TP12  = SQBT  (TE) 

TN12=1./TP12 
TPB=TP*R 

NOTE  THAT  NE(1)  IS  FCR  HCl 

NOTE  THE  INPUT  DATA  FOR  V,VP,U,UP  IS  THE  FOLLOWING; 

RTIP  = 9 V=RCON(J,1)  0 = BCCN(J,2) 

PTYP=10  V = RCON(J,1)  VP  = RCON(J,2) 

BTYP=11  0 = HC0N(J,1)  UP=RCON(J,2) 

DO  300  J=1,NRT 
XLG2=0.0 
NTYP=BTYP  (J) 

GO  TO  (20  1,202,20  3,20  4,20  5,206, 207, 

♦208,209,210,211) ,NTIP 
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201  XLG2*BC0!I(J  ,1)*BC0!I(J  , 3)  *BXP  <RC0II  (J  ,*)/TP) 

?T  OB  HCL  BT  H ATOMS  RTYP=1 

XLG2=V*CONT3*EXP (COMTa/T) 

GO  TO  250 

202  XLG2=BC0N(J  ,2)-*RC0B«J  , 3)  *EXP  fRCON  (J  ,4)/TP) 

VT  ON  U2  BY  H ATOMS  prYP=2 

XLG2=U*CONT3*EXP  (CONTU/T) 

GO  TO  25C 

203  COMTINOE 
GO  TO  250 

204  CONTINUE 

VT  ON  H2  BY  HCL  HOLECDIES  RTYP=4 
VT  ON  H2  BY  AB  ATOMS  RTYP=4 

VT  ON  H2  BY  H2  MOLECCLES  BTYP=4 

DE=  (BE  (2)  -2*RCON  (J , 2 ) ♦B EX E ( 2)  ) 

DE23=DE**P2i 

XLG2=RCON  (J  ,2)*RCON(J  , 3)  *EXP  (PCON  (J  ,4)  *I)E23*T  1 3) 

XLG2  = U*CONT3*EXP  (CONT 4 * DE2 3*T** - . 3 33) 

GO  TO  250 

205  CONTINUE 

VT  ON  HCL  BY  H2  MOLECULES  RTYP=5 
VT  ON  HCL  BY  AP  ATOMS  RTYP=5 

DE=BE  (1)  - 2*FCON  (J,1)*«IEXE(1) 

DE23=DE**P23 

XLG2=BCON(J  ,1)*(FCON(J  ,3)-»RC0N(J  , 5)  * EXP  (RCON  ( J , 6)  *D£  2 3*T  1 3 ) ) 
CE=BE-2VBEXE 

XLG2  = V*  (CONT3+CONT5*i;XP  (CONT6* DE 23*TP* . 3 3 3)  ) 

GO  TO  250 

206  CONTINUE 
GO  TO  250 

207  CONTINUE 

VT  ON  HCL  BY  HCL  MOLECULES  RTYP=7 
DE= (BE (1) -2*RCON  ( J , 1 ) *B EX E ( 1) ) 

D£23=DE**P23 

ILG2=RCON(J  ,1)*(PCON(J  , 3)  ♦ E X P (RCON  ( J , 4 ) * DE2  3*T  1 3)  ♦ 

♦ PCON(J  ,5)  *EXP  (RCON  (J  ,6)*T13)) 

GO  TO  250 

208  CONTINUE 
GO  TO  250 

209  COITINUE 

VV  ON  H2  BY  HCL  MOLECULES  RTYP=9 
H2(0)^  HCL(V-I)  = H2(0-1)>  HCL  (V) 

DE=ABS  (BE  (2) -2*PCON  (J  , 2)  *WEXE  ( 2) -BE  ( 1 ) ♦•  2*RCON  (J  ,1)*BEXE(1)) 
DE23=DE**P23 

XLG2*RCON(J  ,1)*RCON(J  ,2)*(BCON(J  ,3)+RCON(J  , 5)  *EXP  (PCON  (J  ,6)* 
1DE23*T13) ) 

DE  = ABS (BE (2) ♦20BEXE  (2) -BE  (1) -2VBEXE  (1) ) 

XLG2=V*U* (CCNT3+C0NT5*EXP(C0NT6*DE23*TP**-.333) 

GO  TO  250 

210  CONTINUE 

C If  IN  eCL  BTYP*10 
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BCL(f)*  HCL(?P-1)  * HCL<V-1)»  aCL(fP) 

DE=2*(HC0H(J  ,2)-BC0»(J  ,1))*HEXE(1) 

DB23=DE**P23 

XLG2  = BC0!I(J  ,1)*BC0N(J  ,2)*(8CON(J  , 3)  *TN  1 2*  EXP  (BOON  (J  ,U)*DE)^ 
1BCON{J  ,5)  *TP12*BXP  (- (ECON  (J  , 6)  *DE2  3»  BOON  { J ,7))*T13)) 

GO  TO  250 
211  CONTINOE 

I?  IN  H2  BTYP=11 

H2(U)t  H2(0P-1)  = H2(0-1)*  H2  (UP) 

DE=2*(BCON(J  ,2)-BCOB(J  ,1))*HEXB(2) 

DB23=DE**P23 

XLG2  = BCON(J  ,1)*BCON(J  ,2)*(BCON(J  , 3)  * TN  1 2*  EXP  ( BOON  { J ,a)*D5)^ 
1BCON(J  ,5)  *TP12*EXP  (- (BOON  (J  , 6)  *DE2  3^  BOON  ( J ,7))*T1J)) 

250  BX(J)=XLG2 
C NBITE  (6 ,a00)  J,T13, CE23,XLG2 

400  F0BHAT(1X,I3,E15.8,2X,E15.8,2X,E15. 8) 

300  CONTINOE 
BBTOBN 
END 
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C COHSEBVATIOll  OP  HASS 

V2=H/ (A2*RHC2) 

HRITE(6,50) V2 
50  FOBMATC  V2=  • , 1PE20. 8) 

C CHECK  CONVERGENCE 

IP (ABS (V2P-V2)/V2.GT.ERR) GO  TO  100 
EETORN 

200  HRITE(6, 10)T1,T2,P1,F2,V1,V2,Q,H,A2,CP,R 
10  POHSATC  PLON  SUBROUTINE  DID  NCT  CCN  VERGE  • / ( 1 X , 1 P6  E 1 1 . 3 ) ) 

RETURN  1 
END 

PUNCTION  BAD(T,P) 

C THIS  PUNCTION  IS  NOT  YET  OPERATIONAL 

PAD=110. 

HRITE  (6, 1C) 

10  FORHATC  FUNCTION  RAC  IS  NOT  YET  OPERATIONAL,  IT  RETURNS  U=110.') 
RETURN 
END 
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2U1.0  .03282  1.102E-05 
2*2.0  .03296  1.2021-05 
2*3.2  .03312  1.352E-05 


2%*.5  .03330  1.57M1-05 
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